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__ Hadron correlations in jets, ridges and opposite dijets at all px above 2 GeV/c are 

<~| discussed. Since abundant data are available from RHIC at intermediate pt, a reliable 

■^-^ hadronization scheme at that pT range is necessary in order to relate the semihard par- 

tonic processes to the obscrvablcs. The recombination model is therefore first reviewed 
O for that purpose. Final-state interaction is shown to be imijoita.ut for the Cronin effect, 

^ large B/M ratio and forward production. The effect of semihard partons on the medium 

Ch is then discussed with particular emphasis on the formation of ridge with or without 

trigger. Azimuthal anisotropy can result from ridges without early thermalization. Dy- 
, namical path length distribution is derived for any centrality. Dihadron correlations in 

jets on the same or opposite side are shown to reveal detail properties of trigger and 
antitrigger biases with the inference that tangential jets dominate the dijets accessible 
to observation. 
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Hadron Correlations in Jets and Ridges through Parton Recombination 3 

1. Introduction 

Among the many properties of the dense medium that have been studied at RHIC, 
the nature of jet-medium interaction has become the subject of particular cur- 
rent InterestEEl Jet quenching, proposed as a means to reveal the effect of the 
hot medium produced in heavy-ion collisions on the hard parton traversing that 
medium has been confirmed by experiments^ and has thereby been referred to 
as a piece of strong evidence for the medium being a deconfined plasma of quarks 
and gluonsEEMini By the time of Quark Matter 2006 the frontier topic has moved 
beyond the suppression of single-particle distribution at high and into the cor- 
relation of hadrons on both the near side and the away side of jetsinnil The data 
on dihadron and trihadron correlations are currently analyzed for low and interme- 
diate pt, so the characteristic of hydrodynamical flow is involved in its interplay 
with semihard partons propagating through the medium. The physics issues are 
therefore broadened from the medium effects on jets to include also the effect of 
jets on medium. Theoretical studies of those problems can no longer be restricted 
to perturbative QCD that is reliable only at high pT or to hydrodynamics that is 
relevant only at low pT- In the absence of any theory based on first principles that 
is suitable for intermediate px, phenomenological modeling is thus inevitable. A 
sample of some of the papers published before 2008 are given in Refs. p!5]l31j . 

Hadron correlation at intermediate pT involves essentially every complication 
that can be listed in heavy-ion collisions. First, there is the characteristic of the 
medium created. Then there is the hard or semihard scattering that generates 
partons propagating through it. The interaction of those partons with the medium 
not only results in the degradation of the parton momentum, but also gives rise 
to ridges in association with trigger jets and to broad structure on the away side. 
Those features observed are in the correlations among hadrons, so hadronization 
is an unavoidable subprocess that stands between the partonic subprocess and the 
detected hadrons. Any realistic model must deal with all aspects of the various 
subprocesses involved. Without an accurate description of hadronization, observed 
data on hadron correlation cannot be reliably related to the partonic origin of ridges 
and jet structure. 

It is generally accepted that fragmentation is the hadronization subprocess at 
high pt, as in lepton-initiated processes. At intermediate and lower pT recombi- 
nation or coalescence subprocess (ReCo) in heavy-ion collisions has been found to 
be more relevant EIMEIIIII Despite differences in detail, the three formulations of 
ReCo are physically very similar. In Refs. [331 [33] the descriptions are 3-dimensional 
and treat recombination and fragmentation as independent additive components of 
hadronization. In Ref. |35j the formulation is ID on the basis that acoUinear partons 
have low probability of coalescence, and is simple enough to incorporate fragmenta- 
tion as a component of recombination (of shower partons) so that there is a smooth 
transition from low to high p^- Since the discussions on jet-medium interaction in 
the main part of this review are based largely on the formalism developed in Ref. 
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[35j that emphasizes the role of shower partons at intermediate pT, the background 
of the subject of recombination along that line is first summarized along with an 
outline of how non-trivial recombination functions are determined. Some questions 
raised by critics, concerning such topics as entropy and how partons are turned 
into constituent quarks, are addressed. More importantly, how shower partons are 
determined is discussed. 

Large baryon-to-meson ratio observed in heavy-ion collisions is a signature of 
ReCo, since the physical reason for it to be higher than in fragmentation is the 
same in all three formulations HlMIlll The discussion here that follows the formula- 
tion of the recombination model (RM) by Hwa- Yang should not be taken to imply 
less significance of the other two, but only the limits of the scope of this review. 
Considerable space is given to the topics of the Cronin effect (to correct a pre- 
vailing misconception) and to forward production at low and intermediate px in 
Sec. 3. The large B/M ratio observed at forward production cements the validity of 
recombination so that one can move on to the main topic of jet-medium interaction. 

The two aspects of the jet-medium interaction, namely, the effect of jets on the 
medium and that of the medium on jets, are discussed in Sees. 4 and 6, respectively. 
In between those two sections we insert a section on azimuthal anisotropy because 
semihard jets can affect what is conventionally referred to as elliptic flow at low px 
and also because ridge formation can depend on the trigger azimuth at intermediate 
Pt- Much theoretical attention has been given in the past year to the phenomena of 
ridges on the near side and double hump on the away side of triggers at intermediate 
p^ imilolOur aim here is not to review the various approaches of those studies, but 
to give an overview of what has been accomplished on these topics in the RM. The 
focus is necessary in order to cover a range of problems that depend on a reliable 
description of the hadronization subprocess. This review is complementary to the 
one given recently by Majumder,*^^ which emphasizes the region oi px much higher 
than what is considered here, so that factorized fragmentation can be applied. 

Due to space limitations this review cannot go into the mathematical details of 
either the basic formalism or the specific problems. Adequate referencing is provided 
to guide the interested reader to the original papers where details can be found. The 
discussions will mainly be qualitative, thus rendering an opportunity to describe the 
motivations, assumptions and physical ideas that underlie the model calculations. 
For example, the shower parton is an important ingredient in this approach that 
interpolates between what are soft (thermal-thermal) and hard (shower-shower), 
but we have neither space nor inclination to revisit the precise scheme in which 
the shower-parton distributions are derived from the fragmentation functions. The 
concept of thermal-shower recombination and its application to intermediate-pT 
physics are more important than the numerical details. Similarly, we emphasize 
the role that the ridges play (without triggers) in the inclusive distributions of 
single particles because of the pervasiveness of semihard scattering, the discussion 
of which can only be phenomenological. 
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Attempts are made to distinguish our approach from conceptions and interpre- 
tations that are generally regarded as conventional wisdom. Some examples of what 
is conventional are: (1) Cronin effect is due to initial-state transverse broadening; 
(2) large B/M ratio is anomalous; (3) azimuthal anisotropy is due to asymmetric 
high pressure gradient at early time; (4) recombination implies quark number scal- 
ing (QNS) of V2] (5) dijets probe the medium interior. In each case evidences are 
given to support an alternative interpretation. In (4), it is the other way around: 
QNS confirms recombination but the breaking of QNS does not imply the fail- 
ure of recombination. Other topics are more current, so no standard views have 
been developed yet. Indeed, there exist a wide variety of approaches to jet-medium 
interaction, and what is described here is only one among many possibilities. 

2. Hadronization by Recombination 
2.1. A historical perspective 

In the 70s when inclusive cross sections were beginning to be measured in hadronic 
processes the only theoretical scheme to treat hadronization was fragmentation 
for lepton-initiated processes for which the interaction of quark were known to be 
the basic subprocess responsible for multiparticle production. The same fragmen- 
tation process was applied also to the production of high-p^ particles in hadronic 
collisions. Local parton-hadron duality was also invoked as a way to avoid fo- 
cusing on the issue of hadronization.^'^ In dual parton model where color strings 
are stretched between quarks and diquarks, the fragmentation functions (FFs) are 
attached to the ends of the strings to materialize the partons to hadrons, even if 
one of the ends is a diquark.^^ However, at low px in PP collisions quarks are not 
isolated objects in the parton model since there are gluons and wee-partons at small 
so the justification for the confinement of color flux to a narrow string is less 
cogent than at high p^- A more physically realistic description of hadronization 
seemed wanting. 

The first serious alternative to fragmentation against the prevalent scheme for 
hadronization was the suggestion that pion production at low pT in pp collisions 
can be treated by recombination.'^ The simple equation that describes it is 



where Fqq(xi, X2) is the qq distribution, taken to be the product Fq{xi)Fq{x2) of the 
q and q distributions already known at the time among the parton distributions of 
a proton. The recombination function (RF) {xi^X2tx) contains the momentum 
conserving 5{xi + X2 — x) with a multiplicative factor that is constrained by the 
counting rule developed for quarks in hadrons. That simple treatment of hadroniza- 
tion turned out to produce results that agreed with the existing data very well. 

The next important step in solidifying the treatment of recombination is the 
detailed study of the RF. If RF is circumscribed by the characteristics of the wave 
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function of the hadron formed, then it should be related to the time-reversed process 
of describing the structure of that hadron. In dealing with that relationship it also 
becomes clear that the distinction between partons and constituent quarks must be 
recognized and then bridged — a problem that has puzzled some users of the RM 
even in recent years. Since hadron structure is the basis for RF, it became essential 
to have a description of the constituents of a hadron in a way that interpolates 
between the hadronic scale and the partonic scale. It is in the context of filling that 
need that the concept of valons was proposed.^S 

The origin of the notion of constituent quarks (CQs) is rooted in solving the 
bound-state problem of hadrons. However, in describing the structure of a nucleon 
in deep inelastic scattering the role of CQs seems to be totally absent in the struc- 
ture functions such as i'W2{x,Q^). The two descriptions are not merely due 
to the difference in reference frames, CQs being in the rest frame, the partons in 
a high-momentum frame. Also important is that the bound state is a problem at 
the hadronic scale, i.e. low Q^, while deep inelastic scattering is at high Q^. The 
two aspects of the problem can be connected by the introduction of valons as the 
dressed valence quarks, i.e., each being a valence quark with its cloud of gluons and 
sea quarks which can be resolved only by high-Q^ probes. At low the internal 
structure of a valon cannot be resolved, so a valon becomes what a CQ would be 
in the momentum- fraction variable in an infinite- momentum frame. Thus the valon 
distribution in a hadron is the wave-function squared of the CQs, whose structure 
functions are described by pQCD at high Q^. Note that the usual description of 
Q^-evolution by DGLAP has no prescription within the theory for the boundary 
condition at low Q^. That distribution at low is precisely what the valon dis- 
tribution specifies. In summary, the structure function T^{x^Q^) of a hadron is a 
convolution of the valon distribution G^/h{y) and the structure function J^^{z, Q^) 
of a valoiP^ 

T\x, Q2) ^Y.f <dvG.ih{v)T'' {x/y. Q^) , (2) 

where y is the momentum fraction (not rapidity) of a valon in the hadron h. The 
first description of the properties of Gy/h{y) is given in [53 [58], derived from the 
early data J-^{x,Q^). More recent determination of Gy/h{y) is described in Ref. 
[59j where more modern parton distribution functions have been used.l^ 

Gy/ii{y) is the single- valon inclusive distribution in hadron h, and is the 
appropriate integral of the exclusive distribution, G^/^ (2/1,2/2) for pion and 
Gw/p(yi) 2/2) 2/3) for proton. More specifically, Gi,/^(?/i, 2/2) is the absolute square 
of the pion wave function (i'i(?/i)w2(2/2) Itt) in the infinite-momentum frame. Once 
we have that, it is trivial to get the RF for pion (i.e., by complex conjugation), 
since it is the time-reversed process. Thus for pion and proton, we have 

C-TT/ ^ ^1^2 ^ (Xi X2\ , . 

R{Xi,X2,x)= — ^G^,/^(— ,— , (3) 

' \ X X / 
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T3P( \ 2:iX2a;3 (Xi X2 0:3 \ 

BF[xi,X2,x-i,x) ^ ^ — G^lp[ —,—, — ] , (4) 

X \ X X X ' 

where the factors on the RHS are due to the fact that the RFs are invariant distribu- 
tions defined in the phase space element J^idxij Xi, whereas G^jh are non-invariant 
defined in Ilidj/i, as seen in (j2]). The exclusive distribution Gyj^ contains the mo- 
mentum conserving b^Ylii For pion there is nothing else, but for other hadrons 
the prefactors arc given in Rcfs. [59l[6Tj. 

Having determined the RF, the natural question next is how partons turn into 
valons before recombination in a scattering process. Let us suppose that we can 
calculate the multi-parton distribution F{xi, X2) for a q and q moving in the same 
direction, whether at low or high px- If their momentum vectors are not parallel, 
with relative transverse momentum larger than the inverse hadronic size, then the 
probability of recombination is negligible. Relative longitudinal momentum need 
not be small, since the RF allows for the variation in the momentum fractions, just 
as the partons in a hadron can have various momentum fractions. Now, as the q and 
q move out of the interaction region, they may undergo color mutation by soft gluon 
radiation as well as dress themselves with gluon emission and reabsorption with the 
possibility of creating virtual qq pairs, none of which can be made precise without 
a high probe. The net effect is that given enough time before hadronization the 
quarks convert themselves to valons with essentially the same original momenta, 
assuming that the energy loss in vacuum due to soft gluon radiation is negligible 
(even though color mutation is not negligible). For that reason we may simply 
write F{xi,X2)R{xi,X2,x) as multiplicative factors, as done in ([ij, while treating 
F{xi,X2) as the distribution of partons and R{xi,X2,x) as the RF of valons. The 
detail of this is explained in Ref . [5S] . 

The question of entropy conservation has been raised at times, especially by 
those with experience in nuclear physics. In elementary processes, such as q+q — > tt, 
unlike a nuclear process p + n — > d, the color degree of freedom is important. 
Since a pion is colorless, the q and q that recombine must have opposite color. 
If they do not, they cannot travel in vacuum without dragging a color flux tube 
behind them. The most energy-efficient way for them to evolve is to emit soft 
gluons thereby mutating their color charges until the qq pair becomes colorless and 
recombine. Such soft processes leave behind color degrees of freedom from the qq 
system whose entropy is consequently not conserved. It is therefore pointless to 
pursue the question of entropy conservation in recombination, since the problem is 
uncalculable and puts no constraint on the kinematics of the formation of hadrons. 
Besides, the entropy principle should not be applied locally. A global consideration 
must recognize that the bulk volume is increasing during the hadronization process, 
and thus this compensates any decrease of local entropy density. 

After the extensive discussion given above on the RF, we have come to the point 
of being able to assert that the main issue about recombination is the determina- 
tion of the multi-parton distribution, such as Fqq{xi, X2) in ([Tj), of the quarks that 
recombine. Related to that is the question about the role of gluons which have to 
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hadronize also. By moving the focus to the distributions of partons that hadronize, 
the investigation can then concentrate on the more relevant issues in heavy-ion 
collisions concerning the effect of the nuclear medium. 

2.2. Shower Partons 

At low pt in the forward direction the partons that recombine are closely related 
to the low-Q^ partons in the projectile. It is a subject to be discussed in a following 
section. At intermediate and high px the partons are divided into two types: ther- 
mal (T) and shower (S). The former contains the medium effect; the latter is due 
to semihard and hard scattered partons. The consideration of shower partons is a 
unique feature of our approach to recombination, which is empowered by the possi- 
bility to include fragmentation process as SS or SSS recombination. The jet-medium 
interaction is taken into account at the hadronization stage by TS recombination, 
although at an earlier stage the energy loss of the partons before emerging from 
the medium is another effect of the interaction that is, of course, also important. A 
quantitative theoretical study of that energy loss in realistic heavy-ion collisions at 
fixed centrality cannot be carried out and compared with data without a reliable 
description of hadronization. At intermediate pt there is no evidence that fragmen- 
tation is applicable because the baryon/meson ratio would be too small, as we shall 
describe in the next section. 

The fragmentation function (FF), D{x), is a phenomenological quantity whose 
evolution is calculable in pQCD; however, at some low before evolution the 
distribution in x is parametrized by fitting the data. With that reality in mind it 
is reasonable to consider an alternative way of treating the FF, one that builds 
in more dynamical content by regarding fragmentation as a recombination process. 
That is, if we replace the LHS of Eq. ([T]) by the invariant function xDf{x), then the 
corresponding two-parton distribution in the integrand on the RHS is the product 
distributions of two shower partons in a jet initiated by a parton of type i. To be 
specific, consider, for example, the fragmentation of gluon to pion 

xD-^ix) = / ^^SI{x,)Sl (j^) RU^i.X2,x) , (5) 

J X\ X2 \ i X\ J 

where the q and q distributions in a shower initiated by the gluon are the same, but 
their momentum-fraction dependencies are such that if one (xi) is a leading quark, 
the other {x2) has to be from the remainder (1 — xi) of the parton pool. With 
xDg{x) being a phenomenological input, it is possible to solve (jsj) numerically to 
obtain S^{z). It has been shown in Ref. 62J that there are enough FFs known from 
analyzing leptonic processes to render feasible the determination of various shower 
parton distributions (SPDs), which are denoted collectively by Sf with i = q,q,g 
and j — q, s, q, s, where q can be either u or d. If in i the initiating hard parton is 
an s quark, it is treated as q. That is not the case if s is in the produced shower. 
The parameterization of Sf has the form 

Sf{z)^Az''{l~zf{l + cz''), (6) 
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where the dependence of A, a, etc., on i and j are given in a Table in Ref. ^62] , 
Those parameters were determined from fitting the FFs at Q = 10 GeV, and have 
been used for all hadronization processes without further consideration of their 
dependence on pT- It should be recognized that those shower partons are not to be 
identified with the ones due to gluon radiation at very high virtuality calculable in 
pQCD, which is not applicable for the description of hadronization at low virtuality. 

To sum up, in the conventional approach the FF is treated as a black box with a 
parton going in and a hadron going out, whereas in the KM we open up the black box 
and treat the outgoing hadron as the product of recombination of shower partons, 
whose distributions are to be determined from the FFs. Once the SPDs are known, 
one can then consider the possibility that a shower parton may recombine with a 
thermal parton in the vicinity of a jet and thus give a more complete description of 
hadronization at intermediate-pr region, especially in the case of nuclear collisions. 

The SPDs parametrized by (|6]), being derived from the meson FFs, open up 
the question of what happens if, instead of jj' M, three quarks in the shower 
recombine, e.g., uud — > p, or uds — > A. A self-consistent scheme of hadronization 
would have to demand that the formation of baryons is a possibility in a fragmen- 
tation process and that the SPDs already determined should give an unambiguous 
prediction of what baryon FFs are. The calculation has been carried out in Ref. (53] 
where the results for g p and g A in gluon jets are in good agreement with 
dataP^ without the use of any adjustable parameters. To be able to relate meson 
and baryon FFs is an attribute of our formalism for hadronization that has not 
been achieved in other theoretical approaches, and provides further evidence that 
the SPDs are reliable for use at the hadronization scale. 

2.3. Parton distributions before recombination 

In the study of shower partons in a jet we have assumed the validity of the approx- 
imation that the fragmentation process is essentially one dimensional. One may 
question whether the recombination process in a nuclear collision for a hadron pro- 
duced at high pt may necessitate a 3D consideration, since two different length 
scales seem to be involved, one being that of the hadron produced, the other being 
the size of nuclear medium at hadronization time. Indeed, recombination schemes 
formulated in 3D have been proposed, and various groups have independently found 
satisfactory results that are similar to one another IMMMl 

The essence of recombination is, however, not in the 2D transverse plane normal 
to the direction of hadron momentum because if the coalescing parton momenta are 
not roughly parallel, then the relative momentum would have a large component in 
that transverse plane. If that component is larger than the inverse of the hadron size, 
then the two (or three) partons cannot recombine. Thus partons from regions of the 
nuclear medium that are far apart cannot form a hadron, rendering the concern over 
different length scales in the problem inessential. Only coUinear partons emanating 
from the same region of the dense medium can recombine. For that reason the ID 
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formulation of recombination is adequate, as simple as expressed in ([l]). If one asks 
why the relative momentum can be large in the hadron direction, but not transverse 
to it, the answer lies in the foundation of the parton model where the momentum 
fraction can vary from to 1, while the transverse parton momentum kx is limited 
to ~ (hadron radius)"^. The RFs in ^ and Q are related to the ID wave function 
in that framework. 

Having justified the ID formulation of recombination, let us now focus on the 
distributions of the recombining partons at low px, and later at high px- Since 
pQCD cannot be applied to multiparticle production at low px, our consideration 
of the problem is based on Feynman's parton model, which was originally proposed 
for hadron production at low px-^^ In a,pp collision there are valence and sea quarks 
and gluons whose x-distributions at low are known!^ Without hard scattering 
their momenta carry them forward, and they must hadronize in the fragmentation 
region of the initial proton. RM has provided a quantitative treatment of the single- 
hadron inclusive distribution in xp not only for pp, but also for all realistic hadronic 
collisions What is to be remarked here is how gluons hadronize. In an in- 
cident proton as in other hadrons, the gluons carry about half the momentum of 
the host. Since gluons cannot hadronize by themselves, but can virtually turn to qq 
pairs in the sea, we require that all gluons be converted to the sea quarks (thus sat- 
urating the sea) before recombination. This idea was originally suggested by Duke 
and Taylor,!^ and was implemented quantitatively in the valon model in Refs. |58l 
I5T| . The simplest way to achieve that is to increase the normalization of the qq sea 
quarks without changing their x distributions so that the total momentum of the 
valence (unchanged) and sea quarks (enhanced) exhausts the initial momentum of 
the hadron without any left over for gluons. With the Fqq{xi,X2) thus obtained, 
the use of ([T]) results in an inclusive tt distribution that agrees with data in both 
normalization and x spectrumHHlIl Using the appropriate valon distributions of 
pion and kaon, the success extends beyond p tt* to 7r+ — + tt^, — )■ tt*, and 
TT"*" — > in hadronic collisions at low px- In nuclear collisions there is the ad- 
ditional complication arising from momentum degradation when partons traverse 
nuclear medium. It is a subject that will be brought up in Sec. 3.4. 

When Px is not small, then there has to be a semihard or hard scattering at the 
partonic level so that a parton with kx > 5 GeV/c has to be created. In that case 
shower partons are developed in addition to the thermal partons, so the partons 
before recombination can be separated into the following types: TT-|-TS-t-SS for 
mesons and TTT-hTTS-hTSS-hSSS for baryons.^ The thermal partons have kx 
mainly < 2 GeV/c. If one has a reliable scheme to calculate the thermal partons, 
then their kx distributions can, of course, be used in the recombination equation. It 
does not mean that hydrodynamics is a necessary input in the RM. In pA collisions, 
for instance, hydrodynamics is not reliable, yet the Cronin effect can be understood 
in the RM for both proton and pion production without associating the effect with 
initial-state scattering — a departure from the conventional thinking that will be 
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discussed in the next section. In most applications reviewed here, the distributions 
of thermal partons are determined from fitting the data at low px, and are then used 
in the RM to describe the behavior of hadrons at ^ 3 GeV/c. When we consider 
correlation at a later section, careful attention will be given to the enhancement of 
thermal partons due to the energy loss of a semihard or hard parton passing through 
the nuclear medium. It is only in the framework of a reliable hadronization scheme 
can one learn from the detected hadrons the nature of jet-medium interaction, as 
aspired in jet tomography. 

3. Large Baryon/Meson Ratios 

3.1. Intermediate px in heavy-ion collisions 

A well-known signature of the RM is that the baryon/meson (B/M) ratio is large — 
larger than what is customarily expected in fragmentation. The p/n ratio of the FFs, 
i.e., Dp^g(x)/ D^^ij(x), is at most 0.2 at a; ~ 0.3, and is much lower at other values 
of a;.'^ However, for inclusive distributions in heavy-ion collisions at RHIC the ratio 
Rp/TT is as large as ~ 1 at ~ 3 GeV/c, as shown in Fig. 1. Thus hadronization 
at intermediate px cannot be due to parton fragmentation. Three groups (TAM, 
Duke and Oregon) have studied the problem in the Recombination/ Coalescence 
(ReCo) modeElMElMl and found large Rp/.^ in agreement with the data. The 
underlying reason that is common in all versions of ReCo is that for p and tt at 
the same px the three quarks that form the p has average momentum pt/3, while 
the q and q that form the tt has pt/2. Since parton distributions are suppressed 
severely at increasing kx, there are more quarks at px/3 than at pt/2, so the 
formation of proton is not at a disadvantage compared to that of a pion despite 
the difference in the RFs. For either hadron the recombination process is at an 
advantage over fragmentation because of the addivity of momenta. Fragmentation 
suffers from two penalties: first, the initiating parton must have a momentum higher 
than px, and second, the FFs are suppressed at any momentum fraction, more for 
proton than for pion. Thus the yield from parton fragmentation is lower compared 
to that from parton recombination at intermediate px, even apart from the issue 
of B/M ratio. When faced with the question why baryon production is so efficient, 
the proponents of pion fragmentation regard it as an anomaly. Despite efforts to 
explain the enhancement in terms of baryon junction jHEH! the program has not been 
successful in establishing it as a viable mechanism for the formation of baryons.'^ 
From the point of view of ReCo there is nothing anomalous. 

A simple way to understand the px dependence of Rpf-nipr) is to consider the 
ID formulation of ReCo given in Ref. [35J, where the invariant distributions of 
meson and baryon production are expressed as 
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Fig. 1. Com parison of baryon/meson ratio from STAR data to theoretical curves 
ReCoIMEIEl Left panel is from Ref. |66| : right panel is from Ref. | 67| . 



from 



P""^ = / (n|^j ^3,(gi,'Z2,93)i?^(gi,g2,93,p) (8) 

in which all quarks are coUinear with the hadron momentum p. We assume that 
the rapidity y is « 0, so the transverse momenta are the only essential variables, 
for which the subscripts T of all momenta are therefore omitted, for brevity. Mass 
effect at low p-p renders the approximation poor and the ID description inadequate. 
However, in order to gain a transparent picture analytically, let us ignore those 
complications and assume provisionally that all hadrons are massless. Then the 
experimental observation of exponential behavior of the distribution of pions at 
low pt, i.e., dN'^ /pdp oc exp(— p/T), implies that the thermal partons behave as 

where C has dimension (GeV)~^, and given in ^ is dimensionless. When 
thermal partons dominate Fqq and F^q , the multiparton distributions can be written 
as products: T{qi)T{q2) and T {qi)T {q2)T (q^) , respectively. It is then clear from the 
dimensionlessness of the quantities in ([t]) and Q that with the proton distribution 
having dependence, as opposed to the pion distribution being cx C^, the p/ir 
ratio has the property 

dNP/pdp 
^^/-^P^^dN^p'^^P ' 

so long as thermal recombination dominates. This linear rise with p is the behavior 
seen in Fig. 1, although the mass effect of proton makes it less trivial in px- Nev- 
ertheless, this simple feature is embodied in the more detailed computation until 
shower partons become important for > 3 GeY/cW^ 

From the above analysis which should apply to any baryon and meson, it follows 
that the ratios A/K and U/(f> should also increase with pr in a way similar to p/ir. 
Such behaviors have indeed been observed by STAR|IliI2l as have been obtained 
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in theoretical calculation!^ Taken altogether, it means that without TS and TTS 
recombination the B/M ratios would continue to rise with pT- But the data all 
show that the ratios peak at around pt ~ 3 GeV/c. In the RM the bend-over 
is due to the increase of the TS component of the meson earlier than the TTS 
component of the baryon, since two thermal partons in the latter have more weight 
than the single thermal parton in the former. The shower parton distribution S{q) 
in heavy-ion collisions is a convolution of the hard parton distribution fi{k) and 
the S distribution derived from FF, discussed in Sec. 2.2, i.e., 

fi(k) is the transverse-momentum distribution of hard parton i at midrapidity and 
contains the shadowing effect of the parton distribution in nuclear collisions. A 
simple parametrization of it is given in Ref. 74J as follows 

= K- , (12) 

where K = 2.5 and A, B, n are tabulated for each parton type i for nuclear collisions 
at RHIC and LHC. The parameter C, = 0.07 is the average suppression factor that 
can be related to the nuclear modification factor Raa, and was denoted by ^ in 
Ref. |35j and other references thereafter. Since fi{k) has a power-law dependence 
on fc, so does S{q2) on q2 in contrast to the exponential behavior of the thermal 
partons, T{qi). This upward bending of 5(^2) relative to T(gi) is the beginning of 
the dominance of TS and TTS components over TT and TTT components, resulting 
in a peak in the B/M ratio at around pr ^ S GeV/c. Detailed descriptions of these 
calculations are given in Ref. |35J |7S] . 

We add here that the effort made to consider the shower partons before recom- 
bination is motivated by our concern that a hard parton with high virtuality cannot 
hadronize by coalescing with a soft parton with low virtuality. The introduction of 
shower partons is our way to bring the effects of hard scattering to the hadronization 
scale. At the same time the formalism does not exclude fragmentation by a hard 
parton, since SS and SSS recombination at high pT are equivalent to fragmentation 
but in a language that has dynamical content at the hadronization scale. 

One could ask how the RM can be applied reliably in the intermediate-pT region 
before the shower partons were introduced. The approach adopted in Ref. |32j 
does not involve the determination of the hard parton distribution by perturbative 
calculation, but uses the pion data as input to extract the parton distribution at 
the hadronization scale at all pt in the framework of the RM. It is on the basis 
of the extracted parton distribution (which must in hindsight contain the shower 
partons) that the proton inclusive distribution is calculated. Thus the procedure 
is self-consistent. The result is that the p/n ratio is large a,t px ^ 3 GeV/c in 
agreement with data; furthermore, it was a prediction that the ratio would decrease 
as Pt increases beyond 3 GeV/c, as confirmed later by data shown in Fig. 1. 
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3.2. Cronin effect 

The conventional explanation of the Cronin effect i.e., the enhancement of hadron 
spectra at intermediate pt in pA collisions with increasing nuclear size, is that it 
is due to multiple scattering of projectile partons as they propagate through the 
target nucleus, thus acquiring transverse momenta, and that a moderately large-fcr 
parton hadronizes by fragmentation.'^ The emphasis has been on the transverse 
broadening of the parton in the initial-state interaction (ISI) and not on the final- 
state interaction (FSI). In fact, the Cronin effect has become synonymous to ISI 
effect in certain circles. However, that line of interpretation ignores another part 
of the original discoverjl^ where the A dependence of hadrons produced in pA 
collisions, when parameterized as 

— {pA hX) cx , (13) 

dpT 

has the property that apipr) > Q;;t(pt) for all pT measured. That experimen- 
tal result alone is sufficient to invalidate the application of fragmentation to the 



hadronization process, since if the A dependence in ( 13 1 arises mainly from the ISI, 
where the multiply-scattered parton picks up it k^, then the transverse broadening 
of that parton should have no knowledge of whether the parton would hadronize 
into a proton or a pion, so a/j should be independent of the hadron type h. 

A modern version of the Cronin effect is given in terms of the central-to- 
peripheral nuclear modification factor for dAu collisions at midrapidity 

jl/Ng,,) dN'^/prdprjC) 

{llN^^a)dN-/prdpr{P) ' ^''^ 

where C and P denote central and peripheral, respectively, and Ncoii is the average 
number of inelastic NN collisions. If hadronization is by fragmentation, which is a 
factorizable subprocess, the FFs for any given h should cancel in the ratio of ( |l4| ), so 
Rqp should be independent of h. However, the data show that R^p{pt) > R(jp{pt) 
for all PT > 1 GeV/c when C = 0-20 % and P = 60-90 % centralities.^^ See Fig. 
2. Clearly ISI is not able to explain this phenomenon, which strongly suggests 
the medium-dependence of hadronization. The data further indicate that the px 
dependence of Rqp {pt) peaks at pT '^-^ 3 GeV/c for both p and tt, reminiscent of 
the p/tt ratio at fixed centrality in AuAu collisions although the C /P ratio for dAu 
collisions is distinctly different. 

Hadron production at intermediate pt and 77 '--^ in dAu collisions can be 
treated in the RM in a similar way as for AuAu collisions. Although no hot and 
dense medium is produced in a dAu collision, so thermal partons are not gener- 
ated in the same sense as in AuAu collisions, nevertheless soft partons are present 
to give rise to the law-pT hadrons. For notational uniformity we continue to refer 
to them as thermal partons. We apply the same formalism developed in Ref. [3S] 
to the dAu problem and consider the TT+TS-I-SS contributions to tt production 
(TTT-HTTS-fTSS-fSSS for p). The thermal T distribution is determined by fitting 
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the pt spectra at pt < 1 GeV/c for each centrahty; the shower-parton S dis- 
tribution is calculated as before but without nuclear suppression. Unlike the dense 
thermal system created in AuAu collisions, the T distribution in this case is weaker; 
its parameters C and T (inverse slope) that correspond to the ones in Eq. ^ are 
smaller. Furthermore, C decreases with increasing peripherality, while T remains 
unchanged at 0.21 GeV/c.^* Thus thermal-shower recombination becomes impor- 
tant at Pt ^ 1 GeV/c, which is earlier than in AuAu collisions. As a consequence, 
Rqp{pt) becomes > 1 at > 1 GeV/c. That is the Cronin effect, but not due to 
ISI. The same situation occurs for proton production, only stronger.l^ The calcu- 
lated results for the inclusive distributions of both tt and p agree well with data at 
all centralities, hence also Rqp{pt) and R^p{pT)- Fig. 2 shows R^p{pt) for C = 
0-20 % and P = 60-90 % in dAu collisions for h — tt and p;^^ the lines are the 
results obtained in the RMlMZHlThe reason for i?^p > R^p can again be traced to 
3-quark recombination for p and only 2 quarks for tt. When px is large, fragmen- 
tation dominates (i.e. SS and SSS), and both Rcp approach 1, since FFs cancel 
and the yields are normalized by N^oii- No exotic mechanism need be invoked to 
explain the p production process. FSI alone is sufficient to provide the underlying 
physics for the Cronin effect. 
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Fig. 2. Central-to-peripheral ratios for the production of pion (left panel) and proton (right panel) 
in dAu collisions. Data are from Ref. 1771 and lines are from Refs. 1781 179| . 



3.3. Forward production in dAu collisions 

Hadron production at forward rapidities in dAu collisions was regarded as a fertile 
ground for exposing the physics of ISI, especially saturation physicsjMEl since the 
nuclear effect in the deuteron fragmentation region was thought to cause minimal 
FSI. It was further thought that the difference in nuclear media for the Au side 
(77 < 0) and the d side {r} > 0) would lead to backward-forward asymmetry in 
particle yield in such a way as to reveal a transition in basic physics from multiple 
scattering in ISI for 77 ^ to gluon saturation for rj > 0. The observation by 
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BRAHMS^ that Rcp decreases with increasing rj was regarded as an indication 
supporting that viewl^That hne of thinking, however, assumes that FSI is invariant 
under changes in rj so that any dependence on 77 observed is a direct signal from ISI. 
Such an assumption is inconsistent with the result of a study of forward production 
in dAu collisions in the RM, where both Rcp{pt,ii) ^'^d Rb/f{pt) arc shown to 
be well reproduced by considering FSI only.^* Any inference on ISI from the data 
must first perform a subtraction of the effect of FSI, and just as in the case of the 
Cronin effect there is essentially nothing left after the subtraction. 

In Sec. 3.2 the Cronin effect at midrapidity {-q w 0) is considered. The extension 
to 77 > along the same line involves no new physics. However, it is necessary 
to determine the 77 dependencies of the soft and hard parton spectra at various 
centralities. For the soft partons, use is made of the data on dN/drj to modify 
the normalization of T(q, rj) already determined at ry = 0. For the hard partons, 
modified parametrizations of their distributions fi(kT, rj) are obtained from leading 
order minijet calculations using the CTEQ5 pdfl^and the EKS98 shadowing .SS- A 
notable feature of the result is that fi{kT,rf) falls rapidly with kx as rj increases, 
especially near the kinematical boundary kx = 8.13 GeV/c and 7/ = 3.2. Thus TS 
and SS components are negligible compared to TT at large 7/ for any px and any 
centrality, even though the TT component is exponentially suppressed. In central 
collisions there is the additional suppression due to momentum degraduation of the 
forward partons going through the nuclear medium of the target Au. Putting the 
various features together leads to the ratio Rcp{pt,v) shown in Fig. 3(a), where 
the data are from Ref. [52] and the curves from the calculation in Ref. [SJ. It is 
evident that the decrease of Rc'p{pt, v) at > 2 GeV/c as rj increased from 7/ = 
to 77 = 3.2 is well reproduced in the RM. Only one new parameter is introduced to 
describe the centrality and rj dependence of the inverse slope T of the soft partons, 
but no new physics has been added. The suppression of RcpiPTiV) at ?/ > 1 is 
due mainly to the reduction of the density of soft partons in the forward direction, 
where hard partons are suppressed. 

Extending the consideration to the backward region and using the same T{r]) 
extrapolated to 7; < 0, the backward/forward ratio of the yield can be calculated.!^ 
For r] = ±0.75 corresponding to the data of STAR at 0.5 < \r]\ < 1.0 and 0-20% 
centrality,'^ the calculated result on Rb/f for tt"*" + 7r" + p is shown by the solid 
line in Fig. 3(b). While it agrees with the data very well for px < 2 GeV/c, it is 
noticably lower than the data for all charged particles for px > 2 GeV/c. However, 
more recent data on Rb/f{pt) for 7r+ + 7r~ + p + p^^ shown in the inset of Fig. 
3(b), exhibit excellent agreement with the same theoretical curve that should be 
regarded as a prediction. 

The fact that Rb/f is > 1 for all px measured may be regarded as a proof 
against initial transverse broadening of partons, since forward partons of d have 
more nuclear matter of Au to go through than the backward partons of Au. Thus 
if ISI is responsible for the acquisition of px of the final-state hadrons, then Rb/f 
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should be < 1. The data clearly indicate otherwise. 
3.4. Forward production in AuAu collisions 

Theoretical study of hadron production in the forward direction in heavy-ion colli- 
sion is a difficult problem for several reasons. The parton momentum distribution 
at low and large momentum fraction x in nuclear collisions is hard to determine, 
especially when momentum degradation that accounts for what is called "baryon 
stopping" cannot be ignored. Furthermore, degradation of high-momentum partons 
in the nuclear medium implies the regeneration of soft partons at lower x; that is 
hard to treat also. The use of data as input to constrain unknown parameters is un- 
avoidable; however, existent data have their own limitations. Measurement at fixed 
rj cannot be used to provide information ot\ xp dependence unless pt is known.l^ 
Measurement of both 77 and pt has been limited to charged hadrond^ol i^j^at cannot 
easily be separated into baryons and mesons. For these various reasons forward 
production in AA collisions has not been an active area of theoretical investigation. 
However, there are gross features at large r/ that suggest important physics at play 
and deserve explanation. 

PHOBOS data show that particles are detected at 77' > where rj' is the shifted 
pseudorapidity defined by ?y' = — ybeam ■'^ It is significant because it suggests that 
if {pt) is not too small, it corresponds toxp > 1, where xp = {pT/'mp)e^ . Instead of 
violation of momentum conservation, the interpretations in the RM is that a proton 
can be produced in the xp > 1 region, if three quarks from three different nucleons 
in the projectile nucleus, each with Xi < 1, recombine to form a nucleon with xp — 
'^^Xi > 1.1^ That kinematical region is referred to as transfragmentation region 
(TFR), which is not accessible, if hadronization is by fragmentation. The theoretical 
calculation in the RM involves an unknown parameter, k, which quantifies the 
degree of momentum degradation of low-fcy partons, in the forward direction. For 
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K in a reasonable range, not only can nucleons be produced continuously across the 
xf boundary, but also can p/n ratio attain an amazingly large value)^ 

BRAHMS has determined the px distribution of all charged particles at = 
3 2)90. Por ^p^^ _ ^ GeV/c, the corresponding values oi xp for pion and proton 
are, respectively, 0.4 and 0.54. Taking the preliminary value of the p/p ratio at 0.05 
into account, it is possible to estimate the value of n and then calculate the px 
distribution oi p + p + tt^ + 7r~.'^The p/n ratio was predicted to be ~ 1 at ~ 1 
GeV/c. However, at QM2008 the more recent data on Rp/p was reported to have a 
lower value at 0.02 and on Rp/yr a higher value at ~ 4 at pt ~ 1 GeV/c.'^^ Those 
new data prompted a reexamination of the problem in the RM; with appropriate 
changes in the treatment of degradation, regeneration and transverse momentum, 
the very large p/ir ratio can be understood.!^ 

Since p, p and tt production at large rj depends sensitively on q and q distri- 
butions, which in turn depend strongly on the dynamical process of momentum 
degradation and soft-parton regeneration (the parameterization of which requires 
phenomenological inputs), the procedure in Ref. [SS] is to use Rp/yr and Rp/p as 
input in order to determine k and then calculate the x distributions of the hadrons. 
At fixed r] the x and px distributions are related. It turns out that the result on 
the X distribution leads to a large contribution to the px distribution of Rp/T,{px) 
shown by the dashed line in Fig. 4(a). The additional enhancement shown by the 
solid line arises from the mass dependence of the inverse slopes Th due to flow. 
While the ratio is insensitive to the absolute normalizations of the yields, the 
inclusive distribution of all charged particles is not. In Fig. 4(b) is shown the good 
agreement between the calculated result and the data in both normalization and 
shape with no extra parameters beyond n already fixed. 
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It should be noted that the p/i: ratio, shown in Fig. 4(a), is extremely large at 
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77 = 3.2 and modest px < 2 GeV/c. The underlying physics is clearly the suppres- 
sion of q at medium x and the enhancement of p due to 3q recombination, where 
the (valence) quarks are from three different nucleons in the projectile. No other 
hadronization mechanisms are known to be able to reproduce the data on the large 
Rp/ir at large r/. 

3.5. Recombination of adjacent jets at LHC 

So far we have considered only the physics at RHIC energies and the recombination 
of thermal and shower partons, either between them or among themselves. At RHIC 
high-pT- jets are rare, so the shower partons are from one jet at most in an event. At 
LHC, however, high-pT- jets are copiously produced for pt < 20 GeV/c. When the 
jet density is high, the recombination of shower partons in neighboring jets becomes 
more probable and can make a significant contribution to the spectra of hadrons in 
the 10 < pt < 20 GeV/c range, high by RHIC standard, but intermediate at LHC. 
If that turns out to be true, then a remarkable signature is predicted and is easily 
measurable: the p/ir ratio will be huge, perhaps as high as 20W^ 

If a hard parton of momentum kx is produced, shower partons in its jet with 
momenta qi are limited by the constraint qi < kx, so that the recombination of 
those shower partons can produce a hadron with momentum px not exceeding kx- 
However, if there are two adjacent jets with hard-parton momenta kix and k2x, 
then to form a hadron at px from shower partons in those two jets, neither kix 
nor /c2T need to be larger than px, so the rate of such a process would be higher. 
Furthermore, to form a proton at px the shower parton qi can be lower than those 
for pion formation at the same px, so kix can be even lower. Thus i?p/^ in 2-jet 
recombination can be much higher than the ratio in 1-jet fragmentation. 

The probability for 2-jet recombination, however, also depends on the overlap 
of jet cones, since the coalescing shower partons must be nearly coUinear. That 
overlap decreases with increasing kix, so there is a suppression factor in the SS 
or SSS recombination integral that depends on the widths of the jet cones. Using 
some reasonable estimates on all the factors involved, it is found that Rp/^r can be 
between 5 and 20 in the range 10 < < 20 GeV/c, decreasing with increasing 
px^ Although exact numbers are unreliable, the approximate value of Rp/T^ is 
about 2 orders of magnitude higher than what is expected in the usual scenario of 
fragmentation from single hard partons. 

The origin of the large Rp/^r at LHC discussed above is basically the same as that 
for forward production in AuAu collisions at RHIC. In both cases it is the multi- 
source supply of the recombining partons that enhances the proton production. At 
large px at LHC there are more than one jet going in the same direction; at large 
Pl ai RHIC there are more than one nucleon going in the forward direction. In the 
latter case we already have data supporting our view that Rp/^ should be large as 
shown in Fig. 3(a). It would be surprising that our prediction of large Rpf.^ at LHC 
turns out to be untrue. 
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4. Ridgeology Phenomenology of Ridges 

In the previous section the topics of discussion have been exclusively on the single- 
particle distributions in various regions of phase space. Everywhere it is found that 
the B/M ratio is large when pT is in the intermediate range. We now consider two- 
particle correlations, on which there is a wealth of data as a result of the general 
consensus in both the experimental and theoretical communities that more can 
be learned about the dense medium when one studies the system's effect on (and 
response to) penetrating probes. The strong interaction between energetic partons 
and the medium they traverse, resulting in jet quenching, is the underlying physics 
that can be revealed in the jet tomography program.** To calibrate the medium 
effect theoretically, it is necessary to have a reliable framework in which to do 
calculation from first principles, and that is perturbative QCD. Although many 
studies in pQCD have been carried out to learn about the modification of jets in 
dense medium in various approximation schemes )HSE1 they are mainly concerned 
with the effect of the medium on jets at high px, and the results can only be 
compared with data on single-particle distributions, such as Raa{pt)- The response 
of the medium to the passage of hard partons is not what can be calculated in 
pQCD, since it involves soft physics. That is, however, the physical origin of most 
of the characteristics in the correlation data. An understanding of that response 
is one of the objectives of studying correlations. Without the reliable theory to 
describe correlation, especially at low to intermediate px where abundant data exist, 
it becomes necessary to use phenomenological models to relate various features of 
correlation. When all the features can consistently be explained in the framework 
of a model, then one may feel that a few parameters are a small price to pay for 
the elucidation of the jet-medium interaction. 

On two-particle correlation the most active area in recent years has been the 
use of triggers at intermediate or high px to select a restricted class of events and 
the observation of associated particles at various values of r/ and relative to the 
trigger EEU Among the new features found, the discovery of ridges on the near side 
has stimulated extensive interest and activities .1^^ We review in this section only 
those aspects in which recombination plays an important role, which in turn makes 
inferences on the origin of the ridges. We start with a summary of the experimental 
facts. 

4.1. Experimental features of ridges 

The distribution of particles associated with a trigger at intermediate px exhibits 
a peak at small Arj and A(f) sitting on top of a ridge that has a wide range in Arj, 
where A77 and Ac/) are, respectively, the differences of rj and of the associated 
particle from those of the trigger li^l^°*l A 2D correlation function in (Ary, A0) first 
shown by Putschk^^^ at QM06 is reproduced here in Fig. 5(a). STAR has been 
able to separate the ridge (R) from the peak (J), where J refers to Jet, although 
both are features associated with jets. The structure shown in Fig. 5(a) is for 



April 14, 2009 17:22 WSPC/INSTRUCTION FILE Hwa-1 



Hadron Correlations in Jets and Ridges through Parton Recombination 21 

3 < p^'s < 4 GeV/c and p^''°'' > 2 GeV/c in central AuAu collisions. The ridge 
yield integrated over Ar] and A(j> decreases with decreasing iVpart, until it vanishes 
at the lowest A^part corresponding to pp collisions, so R depends strongly on the 
nuclear medium. That is not the case with J. On the other hand, R is also strongly 
correlated to jet production, since the ridge yield is insensitive to p^'^. Thus the 
ridge is a manifestation of jet-medium interaction. 

Putschke further showecP^ that the ridge yield is exponential in its dependence 
on p^^°'^ and that the slope in the semi-log plot is essentially independent of p^'^. 
That is shown by the solid lines in Fig. 5(b). The inverse slope parameterized by 
T is slightly higher than Tq of the inclusive distribution, also shown in that figure. 
Since the px range in that figure is between 2 and 4 GeV/c, we know from single- 
particle distribution that the shape of the inclusive spectrum is at the transition 
from pure exponential on the low side to power-law behavior on the high side. The 
last data point at p^^°'^ — 4 GeV/c being above the straight line is an indication 
of that. Thus the value Tq of the pure exponential part for the bulk is lower than 
what that straight line suggests. The exponential behavior of R should be taken 
to mean that the particles in the ridge are emitted from a thermal source. Usually 
thermal partons are regarded as begin uncorrelated. In the case of R they are 
all correlated to the semihard parton that initiates the jet. We thus interpret the 
observed characteristics as indicating that the ridge is from a thermal source at T, 
enhanced by the energy lost by the semihard parton transversing the medium at 
To. 




Fig. 5. (a) Jet stru cture from Ref. [101| for charged particles associated with a trigger, (b) Depen- 
dence of ridge yielcP^ on p^''°'^ for various . 

The B/M ratio of particles in the ridge is found to be even higher than the same 
ratio of the inclusive distributions in AuAu collisions at 200 GeV. More specifically, 
{p + p)/{tt+ + Tr-) in R for p^"*^ > 4 GeV/c and 2 < < p^'^ is about 1 at 

Pt — 4, GeV/c.'iSli In contrast, that ratio in J is more than 5 times lower. There is 
indication that the A/K ratio in the ridge is just as large.l^ As discussed in Sec. 3, 
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it is hard to find any way to explain the large B/M ratio outside the framework of 
recombination. Since the exponential behavior in pT implies the hadronization of 
thermal partons, the application of recombination very naturally gives rise to large 
B/M ratio, as we have seen in Sec. 3.1. 

Putting together all the experimental features discussed above, we can construct 
a coherent picture of the dynamical origin of the R and J components of the jet 
structure, although no part of it can be rigorously proved for lack of a calculationally 
effective theory of soft physics. There are several stages of the dynamical process. 

(a) A hard or semihard scattering takes place in the medium resulting in a 

parton directed outward in the transverse plane at midrapidity. Because of 
energy loss to the medium, those originating in the interior are not able 
to transverse the medium as effectively as those created near the surface. 
That leads to trigger bias. 

(b) Whatever the nature of the jet-medium interaction is, the energy lost from 

the semihard parton goes to the enhancement of the thermal energy of the 
partons in the near vicinity of the passing trajectory. Those enhanced ther- 
mal partons are swept by the local collective movement outward whether 
or not the flow can be described by equilibrated hydrodynamics initially. 

(c) Since the initial scattering takes place at \ri\ < 0.7, which is the pseudora- 

pidity range of the trigger acceptance, the shower (S) partons associated 
with the jet are restricted to the same range of 77. However, the enhanced 
thermal partons that interact strongly with the medium can be carried 
by the high-77 initial partons that they encounter on the way out and be 
boosted to higher 77. Thus the distribution of the enhanced thermal partons 
is elongated in Ary, but not in because the expansion of the bulk sys- 
tem is in longitudinal and radial directions, not in the azimuthal direction. 
Consequently, the hadronization of the enhanced thermal partons has the 
shape of a ridge. 

(d) In terms of recombination the ridge is formed by TT and TTT recombina- 

tion, while the peak J is formed largely by TS and TTS (or TSS) recom- 
bination, and possibly also by fragmentation (SS and SSS), depending on 
Pt and centrality. Since the J component involves S, it is restricted to a 
narrow cone in A77 and A(j). 

An initial attempt to incorporate all these properties in the RM was made in 
Ref. [in] before the ridge data were reported in QMOB.fi^ By the time of QM08 
ridgeology has become an intensely studied subject, as evidence by the talks in Ref. 

m- 

4.2. Recombination of enchanced thermal partons 

Although the properties of ridges described in the above subsection are derived 
from events with triggers, it should be recognized that ridges are present with or 
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without triggers. That is because the ridges are induced by semihard scattering 
which can take place whether or not a hadron in a chosen px range is used to select 
a subset of events. Experimentally, it is known that the peak and ridge structure is 
seen in auto-correlation where no triggers are used.li^^l The implication is that the 
ridge hadrons are pervasive and are always present in the single-particle spectra. 

Hard scattering of partons can occur at all virtuality Q^, with increasing prob- 
ability at lower and lower Q^. When the parton fcy is < 3 GeV/c, the rate of such 
semihard scattering can be high, while the time scale involved is low enough (~ 0.1 
fm/c) to be sensitive to the initial spatial configuration of the collision system. 
Thus for noncentral collisions there can be nontrivial (j) dependence, which we shall 
discuss in Sec. 5. Hadron formation does not take place until much later, so it is 
important to bear in mind the two time scales involved in ridgeology. Ridges are the 
hadronization products of enhanced thermal partons at late time, which are stim- 
ulated by semihard parton created at early time. In the absence of a theoretical 
framework to calculate the degree of enhancement due to energy loss, we extract the 
characteristics of the thermal distributions from the data. Although hydrodynamics 
may be a valid description of the collective flow after local thermal equilibrium is 
established, it does not take semihard scattering into consideration and assumes 
fast thermaliztion without firmly grounded justification. If the semihard scattering 
occurs in the interior of the dense medium, the energy of the scattered parton is 
dissipated in the medium and contributes to the thermalization of the bulk (B). 
That process may take some time to complete. If the semihard scattering occurs 
near the surface of the medium, its effect can be detected as J + R in these events 
selected by a trigger with the trigger direction not far from the local flow direction, a 
point to be discussed in more detail later in Sec. 4.4. Inclusive distribution averages 
over all events without triggers, including all manifestation of hard and semihard 
scatterings; hence, it is the sum of B + R + J. Since J is associated with the shower 
partons (S), we identify J with the recombination of TS-fSS for the mesons and 
TTS-FTSS-hSSS for the baryons, leaving TT-^TTT for B + R. Thus the exponential 
behavior of the thermal partons is revealed in the exponential behavior of i? -I- i? in 
Pt, for which we emphasize the inclusion of the ridge contribution to the inclusive 
distribution. 

In noncentral collisions the ridges are not produced uniformly throughout all 
azimuthjli^ so dN/ dpT that averages over all (f) has varying proportions of B and R 
contributions depending on centrality. To be certain that we can get a measure of 
the R contribution independent of (j), we focus on only the most central collisions in 
this and the next subsections. Continuing to use the notation fc-r for the transverse 
momentum of the semihard parton at the point of creation in the medium, qt for 
that at the point of exit from the medium, and pT for the hadron outside, we have 
for thermal partons the distribution given in ^ just before recombination. Our 
first point to stress here is that the inverse slope T in (|9]) includes the effects of 
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both B and R. Putting that expression into ([t]) where one takes 

Fqq{qiT,q2T) =T{(11t)T{12t) (15) 

and being more explicit with the RF for pion in (|3]), i.e., 

r>7r/ N qiT<l2Tj;fqiT , 92r 

R {qiT^i2T,PT) ^ — 2 — H ^ 1 ' 1^ 

Pt \Pt pt J 

one obtaind^ 

dN^^" = —e-P^/^ , (17) 

although in 2004 no one was aware of the existence of ridges. From the dataP^ 
for identified hadrons, one can fit the tt"*" distribution for 0-5% centrality in the 
range 1 < pt < 3 GeV/c and get T = 0.3 GeV/c. This value is slightly lower than 
the one given in Ref. |101j which takes the slope of the inclusive distribution in 
the range 2 < pT < 3.5 GeV/c. Ref. |105j provides data for K and p also, which 
have the same value of T as above for 1 < pr < 3 GeV/c, thus confirming that 
the exponential behaviors of the hadronic spectra can be traced to the common 
thermal distribution in ^ through recombination. At lower pT the spectra for K 
and p deviate from exponential behavior because of mass effect, which can largely 
be taken into account by using Et instead of pt, where 

9 1 /2 

Et{pt) = rriT - mo, niT ^ {pt + ml) , (18) 
mo being the hadron rest mass. Thus we write for all hadrons 

^^^^ = ^,(pT)e-^-(''-)/^% (19) 

where At^{pt) — is a constant for pion, but for proton Ap{pT) = A^p^ / p^ 

where Aq is a numerical factor that arises from the wave functions (valon distri- 
bution) of the proton.i^^ Note that the inverse slope is now denoted by , since 
the dataP^ show dependence on hadron type when the distributions are plotted 
as functions of Et- Furthermore, is found to depend on centrality, which is a 
feature that can be understood in the RM as being due to the non-factorizability 
of the thermal parton distributions of uud at very peripheral collisions where the 
density of thermal partons is low .1^221 For central collisions, — 0.35 GeV. We 
summarize the empirical results for tt and p as follows: 

r; = 0.3GeV, (20) 

T^; 0.35(1 - 0.5c) GeV, (21) 

where c denotes % centrality, e. g., c = 0.1 for 10 %. We shall hereafter use to 
denote the inverse slope in Et for B + R, and for that in Et for B only, i.e., 

^ = B,,{pt) = A,ipT)e-^-^P-y^^. (22) 

PTdpT 
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It is hard to find data that describes the bulk contribution only, since the effect 
of semihard scattering cannot easily be filtered out. Indeed, as pr — > 0, there is 
no operational way without using trigger to distinguish B(pt) from all inclusive. 
For that reason the prefactor Ahipr) in (22 1 is the same as that in (19 1. In events 
with trigger above a threshold momentum, semihard partons with lower momenta 
than that threshold can contribute to R; it becomes a part of the background, 
which is experimentally treated as B. Thus the only meaningful way to isolate R 
quantitatively is by use of correlation, while accepting the difficulty of separating B 
and R outside the momentum ranges where the correlated particles are measured. 
Another way of stating that attitude is to accept the experimental paradigm of 
regarding the mixed events as a measure of the background (hence, by definition, 
the bulk), and treating R as only that associated with a trigger. Our cautionary 
point to make is that such a background can contain untriggered ridges. In practice, 
one can take the difference between (19 1 and (22 1 and identify it as the ridge yield 



PrdpT 



Rhipr) = Ah{pT)e' 



-Et{pt)/K' 



where 



rpfl 



1 



ATh. 



n 



AT, 



T, 



(23) 



(24) 



If ATh ^ Th, then the quantity in the square bracket makes a small correction 
to the exp [—Et{pt) /T'^} behavior, and one can determine from the data. The 
only data available that address the ridge distribution are in Ref. |101j where the 
associated particles are in the range 2 < p^^°'^ < 4 GeV/c, exhibiting an approxi- 
mately exponential behavior. It is shown in Ref. |106j by using the data for trigger 
momentum in the range 4 < p%^^ < 5 GeV/c that with ATh = 45 MeV in (24) 
the ridge distributions can be well fitted. The expression for Rhipr) in (23 1 has 
no explicit dependence on p^'^, as is roughly the case with the data. It does have 
strong dependence on which is px in (23 1. Experimental exploration of the 

lower p^^""^ region would provide further validation that (23) needs. The physics 
basis for that distribution is the recombination of thermal partons given in ^ . 



4.3. Trigger from the ridge 

We have discussed above the observation of ridge in triggered events, but to have 
a trigger from the ridge seems to put the horse behind the cart. There must be a 
phenomenological motivation for that role reversal. 

Let us start with the single-pion inclusive distribution that shows an exponen- 
tial decrease in px followed by a power-law behavior. The boundary between the 
two regions is at ~ 2 GeV/c. We have associated the exponential region to TT 
recombination and the power-law region to TS-I-SS. We have also discussed the 
contribution to T from the enhanced thermal partons arising from the medium 
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response to semihard partons. In order to be able to investigate the TT compo- 
nent better without the interference from the shower contribution so that one can 
examine the B + R components cleanly, it would be desirable to be able to push 
the TS+SS components out of the way. That is not possible with the light u and d 
quarks, but not impossible with the s quark, since the heavier quark is suppressed 
in hard scattering. 

If one observes the hadrons formed from only the s quarks, either </> or f2, one 
finds exponential behavior at all px measured, which in the case extends to 
as high as 5.5 GeV/c.'^^^'l^^ The absence of any indication of up-bending of the 
distributions clearly suggests that the source of the s quarks is thermal in nature 
and that no shower partons participate in the formation of and f2. That problem 
is studied in Ref. |114l along with K and A production. Indeed, the data can be 
well reproduced by TT for TTT for fl, TT+TS for K, and TTT+TTS+TSS for 
A. 

Since s quarks in S make insignificant contribution to f2 production for px < 6 
GeV/c, and since thermal partons are uncorrelated, it is reasonable to expect that 
the fl observed has no correlated particles. It was therefore predictecP^^ that if fl 
in 3 < pt < Q GeV/c is treated as a trigger particle, there should be no associated 
particles appearing as a peak in A0 on the near side. Within a year that prediction 
was falsified by a report at QM06 showing that there is a near-side peak after 
background subtraction in A0 for in the range 2.5 < p^'^ < 4.5 GeV/c in central 
AuAu collisions with charged particles in the range 1.5 < p?^^""^ < p^'^ .EH The data 
created a dilemma: is f2 created by a jet or not? If it is, why is the px distribution 
strictly exponential with no hint of jet characteristics? If it is not, why is there a 
Acj) peak in azimuthal correlation? The dilemma became known as the fl puzzle .l^^ 

The resolution of that puzzle is in the recognition that both the trigger f2 and 
the associated particles are in the ridge, first conjectured in Ref. |117j and later 
quantified in Ref. |118j . Jets are involved, since without jets there can be no ridge. 
But not all jet structures exhibit a prominent peak above a ridge. It depends on the 
trigger particle and the ranges of p^'^ and p^^°'^. Consider the jet yield compared to 
the ridge yield at 3 < p^'*^ < 6 GeV/c in AuAu coUisions at 0-10% centralityEiliiHI 
The J/R ratio at p^^""^ ~ 1.2 GeV/c decreases as the trigger particle goes from h 
to Kg and then to A/A. For A/A trigger and unidentified charged h associated, the 
J/R ratio is < 10~^ for IAt^Ij < 0.7 and \Arj\u < 1.7. Since the A/A trigger must 
contain an s quark which is absent in the shower, the participating s quark must be 
a thermal parton. For p^'® near 3 GeV/c, thermal s quark around 1 GeV/c or less 
can be quite abundant, and the initiating semihard parton need not be very hard. 
With p^^""^ as low as ~ 1 GeV/c, the light hadrons in R dominate over those in J, so 
J/R is small. As the strangeness in the trigger increases, more thermal s quarks are 
involved with less dependence on shower. J/R is likely to be even smaller, although 
present data with trigger lack statistics to show the A77 distribution. If the jet 
structure shows mainly a ridge with negligible peak (J) in A77, we have referred to 
it as a phantom jet,'ii3i.e., a Jet-less jet. The corresponding A(/) distribution should 



April 14, 2009 17:22 WSPC/INSTRUCTION FILE Hwa-1 



Hadron Correlations in Jets and Ridges through Parton Recombination 27 

then be dominantly R. Since the initiating semihard partons are either gluon or hght 
quarks, the usual jet structure may still be seen if the trigger particle is ordinary. 
But for f2 trigger the structure is very different. 

Now, we can address the f2 puzzle. The enhanced thermal partons generated by 
the semihard parton contain s as well as u and d quarks. Three s quarks anywhere 
in the ridge can recombine to form a trigger il. Other quarks, in particular the light 
quarks, in the ridge can form associated particles. The pool of enhanced thermal 
partons are all correlated to the semihard parton in every event selected by the 
trigger, so the associated particles are all restricted to |A(/)| < 1. Since the s quarks 
that form the fl are thermal, the f2 spectrum in px is exponential. Hence the 
puzzle is solved. The trigger can be from the ridge. 

The above description outlines a detailed calculation of the Acf) distribution of 
charged hadrons produced in association with il}^^ The background is calculated 
in the RM using previous parametrizations, and the height agrees with the data. 
That is important, since the ridge signal is less than 4% of the background height. 
Two adjustable parameters are used to fit the A<j) distribution of the ridge, but 
then the yield/trigger is calculated as a function of p^'® without further unknowns 
in the model, and the result is in agreement with the datapnUU as shown in Fig. 
6. The solid and dashed curves in the two panels of that figure are the results of 
the calculation using two values of the strength of enhanced thermal partons that 
differ by only 1%, yet the height of the ridge varies by about 20%. That is because 
the ridge is the difference between large numbers of B + R and B. Such accuracy 
is beyond the scope of any dynamical theory to achieve. The phenomenological 
approach adopted in Ref. |118j has been the only one that offers a quantitative 
understanding of the problem. 




3 4 
p^"9 (GeV/c) 



Fig. 6. (a) Left panel; DataT^on associated particle distribution in A4> for three hyperon triggers 
at 2.5 < py'*^ < 4.5 GeV/c and 1.5 < pf^°'' < p^"*^. (b) Right panel: Near-side ridge yield 
associated with Q trigger. The lines are from calculations in the RK^E^ for particles associated 
with Q trigger, with the solid and dashed lines for slightly different normalization constants of the 
Gaussian peak. 
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Recently, PHENIX has shown datcPSl that can be interpreted as support for 
the notion of trigger from the ridge. At p^'^ < 4 GeV/c (for unidentified trigger) 
the per-trigger yield of the associated particles is found to be less than expected 
from fragmentation, and the "dilution" effect is attributed to the increase of the 
number of triggers due to soft processes. In our language fragmentation is the SS 
component, and the TS component due to medium effect gives an increase already 
over fragmentation for pt < Q GeV/c. The additional dilution effect at < 4 
GeV/c is due to TT recombination, which is enhanced by the ridge contribution. 

4.4. Dependence of ridge formation on trigger azimuth 

So far our consideration in ridgeology has been concerned mainly with the depen- 
dence on px- Now, let us turn to azimuthal correlation, although our discussion of 
the main topic of azimuthal anisotropy of single-particle distributions is deferred 
to the following section. Our present focus is on the correlation between the direc- 
tions of the trigger and ridge particles in the transverse plane, a topic that can be 
discussed prior to anistropy because it is mainly a problem in ridge formation. 
The subject was stimulated by the reporifISS that the ridge yield depends on the 
azimuthal angle (jjg between the trigger angle 4>t and the reaction plane rp, even 
for nearly central collisions, but especially for noncentral collisions. 

Since geometry is an important factor that influences the 4>s dependence, it is 
necessary to treat carefully the initial configuration of the problem: (a) the point of 
origin {xo,yo) of the semihard parton in the almond-shaped overlap region, (b) the 
angle 0s of the parton's trajectory, (c) the density of the medium, D{x,y), along 
that trajectory, and (d) the point of exit from the medium. In the approximation 
that the medium does not expand very much during the time that the semihard 
parton near the surface traverses the medium, it is not difficult to calculate the 
path length, but it is much more diflicult to calculate the energy loss that depends 
on D(a;, y). Even if there is a reliable way to account for the effect of the medium 
on the semihard parton, there is no known way to translate that to the effect of 
the parton on the medium. The enhancement of the thermal partons that lead to 
the ridge particles takes time to develop, during which the medium expands. Local 
flow direction depends on where the enhanced thermal partons are in the overlap, 
which evolves into elliptical geometry. Each of the various parts of the process can 
be represented by a factor that can be expressed in terms of variables that have 
reasonable physical relevance. Without entering into the details that are described 
in Ref. |121j we can outline the essence here. 

Let P(xQ,yQ,t) denote the probability of detecting a parton emerging from the 
medium, where t is the path length measured from the initial point (a;o,?/o) to the 
surface along a straight-line trajectory at angle ^g. Let C{'il){x, y),4>s) be a function 
that describes the correlation between the enhanced thermal partons along the 
flow direction ^p{x,y) and the semihard parton direction 0^. Finally, let T{x,y,(j)) 
describe the fluctuation of the angle cj) oi a. ridge particle from the average flow 
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direction. For fixed {xo-,yo) and (ps, the ridge-particle distribution in (p is then 



Jo 

where A'^ is a normalization constant related to the rate of production of the ridge 
particle and ^ is the fraction of the path length t along the trajectory starting at 
at (xq, yo) and ending at 1. For the observed distribution it is necessary to integrate 
over all (xo,yo)- Not every semihard parton included in that integration gets out 
of the medium to generate a particle that triggers the event. The ridge distribution 
per trigger is therefore normalized by the probability of the ridge-generating parton 
emerging from the medium 

, , s J dxodyoR{(t),(j)s,xo,yo) ,„„x 
J dxodyoP[xo,yo,t) 

where the integration is over the entire region of initial overlap. There is no explicit 



dependence on p^'^ and p^^°'^ in (25 1, but the parameters specifying the different 
factors in the equation do. We leave the momenta fixed in the ranges 3 < p^'^ < 4 
GeV/c and 1.5 < p^^""^ < 2.0 GeV/c, as specified in the experiment and focus 
only on the dependencies on and 0^. 



The most important factor in (25 1 is C(x,?/, 0^), which is parametrized as 

\2' 



C{x,y,(t)s) = exp 



2A 



(27) 



where A specifies the Gaussian width of the correlation between the directions 
of semihard parton (pg and local flow ip{x,y). If ^{x,y) is close to (ps for most 
of the points {x,y) along the trajectory, then the thermal partons enhanced by 
successive soft emissions are carried by the flow along the same direction. The 
effects reinforce one another and lead to the formation of a ridge in a narrow 
cone. On the other hand, if ■0(a:, y) is very different from 0^, then the enhanced soft 
partons are dispersed over a range of surface area, so their hadronization leads to no 
pronounced effect. That is the essence of the correlated emission model (CEM).!^^ 
The probability P(xo,yo,t) depends, in addition to the nuclear thickness at 
(a;0)2/o)i the survival probability S(t), which is exponentially behaved in t. Most of 
the semihard partons that can emerge from the medium are created in a layer near 
the surface. In Ref. |121j the thickness of that layer is set to be approximately Ra/4:. 
A more detailed study of that will be described in Sec. 6 below. The fluctuation 
distribution T{x,y,(j)) turns out not to influence the final result in any sensitive 
way. 

A pictorial description of the origin of (ps dependence of ridge yield in GEM is 
shown in the left panel of Fig. 7, where (a) 0^ ~ and (b) cps ^ 70° in noncentral 
collisions. The trigger directions (j)s are represented by the thin arrows, while the 
flow direction ip(x,y) are depicted by the thick arrows that are normal to the 
surface. When the two types of arrows match, the reinforcement leads to a strong 



April 14, 2009 17:22 WSPC/INSTRUCTION FILE Hwa-1 



30 Rudolph C. Hwa 



ridge. In (a) the matching condition is met where the density is higher than in 
(b) where the matching arrows occur near the top of the overlap. Thus the former 
case for (/i^ ~ has higher ridge yield than in the latter case for larger (pg- That 
difference becomes more drastic as the collision centrality becomes more peripheral. 
The results on the ridge yield integrated over (p are compared to the data in the 
right panel of Fig. 7. In (a) of that panel the normalization is adjusted to fit the data 
point at the lowest (pa and the shape of the (ps dependence is adjusted by varying 
the parameter A in (27l. The value determined is A = 0.11; it corresponds to a 
correlation cone of half- width ^ 20° . The curve in part (b) of that panel for 20-60% 
centrality is obtained without any further adjustment. Although the agreement with 
data is not perfect, the curve does reproduce the trend of steep descend, followed 
by a region of flatness or even a small up-bending. The behavior at large <ps is due 
to the additional contribution of semihard partons from the left-half ellipse at high 
Do- 
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Fig. 7. Left panel; schematic sketches of trigger di rect ions (thin arrows) and flow directions (thick 
arrows) in noncentral collisions. Right panel: Datd^ ^ on the ridge yield/trigger vs 4>s for (a) 0-5% 
and (b) 20-60% centrality. The lines are from CEIVPH with common normalization adjusted to 
fit the left-most data point in (a). 



Without being integrated over <p, the ridge distribution R{(p,(ps) given by (25) 
and (26 1 describes the dependence of the ridge yield on A0 = (p — (pg. An interest- 
ing discovery upon careful study of R{(p,(ps) is that its dependence on is not 
symmetric across the = point, depending on the sign of (ps- In Fig. 8(a) is 
shown its behavior for six positive values of <ps- Not shown are the mirror images 
of those curves across = for corresponding negative values of (ps- Since the 
measurementjiSil averages over both positive and negative values of cpg, there is no 
observable asymmetry in the data. The reason for the asymmetry in Fig. 8(a) is 
that for (pg > Q most of the points along the trajectory have ip{x,y) < <pg, since 
tp{x,y) is generally normal to the surface. Hence, the ridge particles are mostly at 
A(p < 0. The reverse is true for cpg < 0. The left shift of the peaks in A(p in Fig. 
8(a) has been confirmed by STAR recently.t^^Il 
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Fig. 8. (a) Left shift of the ridge yield in A0 for 
function defined in Eq. (j29|. 
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Is > in CEM.E^ (b) Inside-outside asymmetry 



Another way to test the asymmetry is to measure a quantity A((j)s) called inside- 
outside asymmetry function.li^il To that end, define for < 0s < 7r/2 

r+(0,)- d0i?(</),0,), Y^icj),)^ r^^^ dcj)R{c^,c^,) (28) 

and for — 7r/2 < (f>s < reverse the definition. Then for any (jig define 

This asymmetry function should always be positive whether there is a left shift for 
(ps < or a, right shift for 0^ > 0. By reflection symmetry it vanishes at 0^ = and 
±7r/2. In CEM the properties of A{(t)s) at two centralities are shown in Fig. 8(b). 
STAR reported very recently at QM09 that such an asymmetry has indeed been 
found to exist in the data on the ridges for various ^sEHElll 



5. Azimuthal Anisotropy 

The azimuthal dependence of single-particle distribution has been studied ever since 
the beginning of heavy-ion physics.^ ^^ l ^^^ l Hydrodynamical model at low pjCsHUSll 
and ReCo at intermediate pj] 33 | 34 | i33 | i34 | j^^^yg been able to describe the data on 
elliptic flow very well. The only points worthy of comments here are those outside 
the realm of what has been covered in the references given above. There are then 
only three points: (a) Is early thermalization necessary? (b) What is the role of 
the shower partons? (c) At what px does the quark number scaling begin to break 
down? 



5.1. Effects of ridge formation at low px 

The usual hydrodynamical treatment of elliptic flow assumes rapid thermalization 
with initial time of expansion set at tq = 0.6 fm/c. Such an early time of equi- 
libration has never been shown to be the consequence of any dynamical process 
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that is firmly grounded and commonly accepted. The question then is whether the 
azimuthal anisotropy can be driven initially not by pressure gradient at t < 1 fm/c, 
but by some other mechanism that is sensitive to the early spatial configuration. 
That mechanism is suggested in Ref. |135j to be semihard scattering, the rate of 
which can be high for parton fcr around 2-3 GeV/c, while the time scale involved is 
low 0.1 fm/c). Semihard partons created near the surface of the nuclear overlap 
can lead to a continuous range of ridges that is shaped by the initial geometry. The 
effects of such ridges are not considered in the usual studies in hydrodynamics, but 
should not be ignored. It is shown in Refs. |1351I106] that the azimuthal anisotropy 
of the ridges produces the observed V2 at low pT without the use of hydrodynamics. 
Elliptic flow at some point of the expansion of the system may well be describable 
by hydrodynamics, but early thermalization is not required. 

For dN /d<j} of single-particle distribution, no triggers are used and ridges are 
generated by many semihard partons produced in each event. Semihard scatterings 
take place throughout the overlap, but only those occurring near the surface and 
directed normal to the surface can lead to the development of ridges, as discussed 
in Sec. 4.4. If a drastic simplification is made to require ridge particles to be also 
directed only in the the directions normal to the almond-shaped initial boundary 
in the transverse plane, then the distribution of all ridges, R{pt, 0), is restricted 
to the region cf) E TZ, which is a set of angles defined by 

101 < $(6) and Itt - 01 < $(&) , (30) 

where <i>(6) — cos"^(6/2i?^).'i^ Fluctuation from the restricted range is, of course, 
possible, not only because the semihard partons can have any scattered angle, but 
also because the ridge particles can fluctuate from the directions of the parton tra- 
jectories. It is shown in the Appendix of Ref. |106j that the region TZ is nevertheless 
a good approximation even when all those effects plus elliptic geometry are taken 
into account, provided that the inaccuracy in the regions around |0| ~ 7r/2 for non- 
central collisions is unimportant — which is indeed the case, since the density at 
the upper and lower tips of the ellipse is low so ridge production there is suppressed. 
Thus in the box approximation the ridge distribution is 

i?(pr,0) = i?(PT)e(0) , (31) 

where 

6(0) = 0(<i> - 101) + 0(<i> - k - 01) . (32) 

In the assumption that the above anisotropy from the ridge is the only dependence 
in dN/d(j}, replacing the usual assumption of rapid thermalization and the conse- 
quent pressure gradient, the bulk medium is then isotropic and the single-particle 
distribution at low pT can be written in the form 

- BiPT) + i?(PT)e(0) . (33) 

PTdpTUCp 
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The normalized second harmonic in (j) can then be calculated analytically, yielding 

sin2$(6) 

MvtM = (COS20), = ,^(^^)/^(J)V2a>(6) ■ (^4) 



At low pt the first term in the denominator is much larger than the second, so (34) 
is reduced to the even simpler formula 

v^iprM^'^^^.^^^mh) . (35) 
ttB{pt) 

This is such a compact formula that its validity should be checked regardless of its 
derivation for the sake of having available a simple description of what is usually 
called elliptic flow after extensive computation. 



The Pt distributions of B{pt) and R{pt) are given, respectively, in (22 1 and 



(23 1. The latter can be written in a form that factors out B{pt) so that R/B has 



the simple form 



R{pt) 

B{pt) 



e 



which becomes Et{pt) /T'fl at low pt, since T^' is large. It then follows from (35) 
that for small Et 

v^{ET,b) ^ ^sm2m ■ (37) 

Thus the initial slope in Et depends only on T," (which sets the scale) and the 
geometric factor sin2<i>(6). That factorizability is in agreement with the data for 
pion, but not so well for proton for mid-central to peripheral collisions. That is 
because T^' develops a b dependence when the centrality c is above 0.3, as can 



be seen in (21 1. While is essentially constant T^' (c) can be approximated by a 
polynomial in c^^^ 

r;' = 1.7 GeV, T" = 2.37 (1 - 1.05c + 0.260^) GeV . (38) 



For Et ~ 1 GeV the full expression in ( 36 1 should be used in ( 35 1 , instead of the 



small Et approximation in (37). The results for tt and p are shown in Fig. 9, and 
reproduce the dataP^ very well. The reason why Tp(c) decreases with increasing 
peripherality, resulting in similar trend in Tp'(c), is that 3-quark recombination is 
more difficult when the thermal partons become less abundant at lower medium 
density as the collisions get more peripheral. At fixed pt ^ 0.5 GeV/c the b depen- 
dences of ^2 for h = TT and p are studied in Ref. [106] in which it is shown that the 



characteristics of the data are well captured by the simple formula in ( 34 1 



5.2. Effects of shower partons at intermediate px 

As Pt is increased to above 2 GeV/c, it is necessary to consider the role played 
by the shower partons, which introduce <j) dependence due to jet quenching of hard 
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parton that depends on path length in the medium. The dominance of TS + SS 
recombination over TT leads to a change in the px dependence of V2{pt)- 



The shower parton distribution S{q), given in (111, due to a hard parton pro- 
duced at momentum k, is for central collisions and is averaged over all </). Now 
for noncentral collisions with anisotropy, that formula needs to be generalized. 
Assuming that the energy loss of a hard parton is proportional to \/P where k' is 
the initial parton momentum '34' 136 ' ^.^^ write Afc = fc' — fc in the form 



Afc = €{b)i{h, 



(39) 



where £(6, (f) is the normalized path length and e(6) is the energy-loss coefficient 
that depends, apart from geometrical factors, a parameter eg that is to be de- 



termined. After solving (39) for fc', and replacing Cfii^) in (111 by fi{k') at the 



shifted momentum, one can keep the first two non- vanishing terms in the harmonic 
expansion of fi{k'{k, 6, cf))) and get 

/,(fc'(fc, 5, ^)) = h{k) [go(fc, b) + 232 (fc, h) cos 24>] , (40) 

where go and §2 can be determined explicitly in terms of eo, fc and The 
shower contribution to the single-pion distribution, when averaged over all 0, may 
be written in the symbolic form 

d7Vj^+'5'5(6) 1 



, -go{Kh)m®{TS + SS) 

prdpT 27r 

The contribution of the shower component to (pti b) is 

g2{k,b)Mk)(E>{TS + SS) 
go{k,b)f,{k)(g>{TS + SS)- 



W2''^(PT,&) 



(41) 



(42) 



The thermal component V2'^^{pT,b) is as given in (34l and (36). The overall V2 is 
obtained from the above with the help of an interpolating function W{pt, b) 



V^ipT, b) = «2"'"'(PT, b)W{pT. b) + vr\pT, b) [1 - W{pT, b)] 



(43) 
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where 

with TT being the thermal, and TS + SS the shower, components of the (/)-averaged 
(INt^ /pTdpT- By fitting the single-pion distribution at 0-10% centrality over the 
range 1 < < 6 GeV/c using all three TT + TS + SS terms, eq is found to be 0.55 
GeV^/^. It is then possible to calculate V2{pT,b) without any further adjustment; 
the result is shown in Fig. 10(a). The data points for Et > 1 GeV are from Ref. 
[137j . The saturation of in that range is thus interpreted in the RM as being due 
to the shower partons, where W{pT,b) is suppressed and 52 is much smaller than 
gQ at high k. 




For proton V2 the general procedure in the calculation is similar to that for 
pion, except for an additional thermal parton to incorporate. There is also the 
complication of b dependence in the inverse slope in (21 ). Taking them all into 
account the result for V2 is shown in Fig. 10(b). Comparison with dataEMUZl is 
acceptable, although more accurate data are needed to check the calculated results 
at high Et and b. 

What is learned from this study is that the main source of ((> anisotropy is the 
path-length dependence of jet quenching, which is parametrized by one unknown 
eo that is determined by fitting the inclusive distribution at one value of pt- The 
characteristics of V2{pT,b) in Fig. 10 are obtained without any more adjustable 
parameters. 

It should be noted that we have refrained from using the term "elliptic flow", 
except in reference to past work based on hydrodynamics. In Sec. 5.1 the empha- 
sis is on the effect of ridges due to semihard scattering that are not taken into 
account by early-time hydrodynamics, although the exponential behavior of the 
thermal partons may well be the result of late-time hydrodynamical flow. Then in 
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this section hard scattering is incorporated in the treatment of jet- medium interac- 
tion through TS recombination before fragmentation dominates. Such interaction 
is outside hydrodynamics, so the overall characteristics of V2{pT,b) are not the 
properties of flow. In Ref. [M] the effects of jet quenching on the hard partons are 
also considered, but since only the fragmentation of those partons is included, the 
transitional contribution from TS interaction is not explicitly taken into account. 
That turns out to be important in the intermediate pt region. 



5.3. Breaking of quark number scaling 

Quark number scaling (QNS) has long been regarded as a signature of 
recombinatiorEmmill and has been verified in a number of experiments ESlElIElZl 
It has been regarded as a statement of the universality oi V2{pt / nq) / Uq, where Uq is 
the number of constituent quarks in the hadron h. It is based on the assumption of 
factorizability of the distribution of the quarks that recombine, i.e., the multiquark 
distributions in ([7| and ^ (but with (pi dependencies included) can be written as 

rig riq 

i^„, (gi,,/)i;g2,<^2;---) ^n^.fe,'/'.) = n^»fe) [l + 2v^{q,) cos icj,,] . (45) 

i=l 1=1 

Coupled with the assumption that the RF has the simple form d{qi — px/riq), it 
then follows trivially from ([t]) and (|8| that 

«2(pt) = nqV^ipr/Tiq) , (46) 



if the V2{qi) of all quarks in (45 1 are the same, denoted by V2{q)- 

From our discussions throughout this paper it is clear that none of the 
above assumptions are valid under close examination. TT+TS+SS for pion and 
TTT-hTTS-hTSS-hSSS for proton are obviously not factorizable. Even at low pr 
where only the recombination of thermal partons is important, the inverse slopes 



r^, given in (20 1 and (21 1, are not the same for tt and p. Consequently, the R/B 



ratios for tt and p are different, as seen in (36) and (38 1, resulting in different 
V2{Ej',b). Furthermore, the wave functions of tt and p are very different, since the 
pion is a tightly bounds state of the constituent quarks, while the proton is loosely 
bound. That means the momentum fractions of the quarks (valons) are not 1 /2 for 
pion and 1/3 for proton. It is then a very rough approximation to write the momen- 
tum conservation (5-function, SCSqi — px), as S^Uqqx — Pt) with a common qx- At 
intermediate px where shower partons become important, we have seen that they 



acquire the (j) dependence of the hard parton, given in (40 1, so (^i) for the shower 
is different from V2{qj) for the thermal parton. Even if TS and TTS contributions 
dominate, one can at best, by ignoring all other complications, have 

(pt) = v^iqi) + vi{q2), v^{pt) = vl{qi) + «J(g2) + v^iq^) . (47) 



They would not lead to QNS, as expressed in (46 1. 

Most data in support of QNS are for minimum bias and at low Et- It is shown 
in Ref. [1061 that the calculated result at Ex/nq < 0.5 GeV does exhibit QNS in 
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agreement with the data, but the scahng is broken above that. The breaking of 
QNS is due primarily to the nonequivalence of the ((> dependencies of the thermal 
and shower partons. The important point to stress here is that the breaking of 
QNS at intermediate Et does not imply the failure of recombination (in fact, it is 
expected), but the validity of QNS at lower Et does confirm recombination as the 
proper hadronization process. 



6. Hadron Correlation in Dijet Production 

In Sec. 4 the subject of discussion is ridgeology, which is the study of the effects of 
jets on the medium. In this section we consider the reverse, i.e., the effects of the 
medium on the jets for px not so extremely high as to exclude TS recombination. 
In triggered events the ridge is the broad pedestal on top of which sits the peak. 
The structure of that peak, when compared to the jet structure in pp collisions, 
reveals the medium effect on hard partons. Since the peak is restricted to a narrow 
cone around the trigger direction, it is natural in the RM to associate it with 
shower partons, while the ridge being broad is associated with the thermal partons. 
Hadrons closely correlated to a trigger particle as jets exhibit peaks on both the 
same and away sides .^^^11^^ How the structures of the jets on the two opposite sides 
differ from each other is a strong indication of the difference in energy losses in the 
two jets, since their path lengths in the medium are generally different. In realistic 
collisions even when centrality is chosen to be within a narrow range, the path 
lengths can vary significantly depending on the location and angle of a scattered 
hard parton. Thus a careful study of the properties of the near- and away-side jets 
must start with finding a good description of the variation of energy loss within 
each class of centrality. 



6.1. Distribution of dynamical path length 

By dynamical path length we mean not only the geometrical length of a trajectory, 
but also the medium effect along that trajectory. To quantify that in an analytic 
expression, we need to revisit the single-particle distribution discussed earlier for 
central collision, but now formulated in a way appropriate for any centrality. For 
pion production we have for y ^ 



pdp 



C2 

6 



■,-p/T 



p2 L-j 



dq 



p 



q 



(48) 



where all momenta are in the transverse plane with the subscript T omitted. The 



first term on the right side is from (17) for TT recombination; the centrality de- 
pendence of C is given in Ref. [143J. The first term in the square bracket is for TS 
recombination which will be detailed below, and the second term is the fragmen- 
tation function that is equivalent to SS recombination. Fi{q) is the distribution of 
parton i at the surface of the medium with q denoting the momentum of the hard 



parton there. It differs from the distribution fi{k) given in (12 1, that describes the 
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hard parton with momentum k at the point of creation. In Sees. 2 and 3 where only 
central collisions are considered, an average suppression factor Q is used, as shown 



F^{D^ / ? / dkkMk)G{q,k,t), (49) 



in (11 ). We now replace C/i(^) by a path-dependent term related to Fi(q) by 

rL 

L 

where G{q, k, t) is a degradation factor that describes the decrease of parton momen- 
tum from fc to g as the hard parton traverses a distance t through the medium. L 
is the average maximum length of that trajectory. In the limit L — s- 0, Fi{q) should 
become the parton distribution function Fi(k) for pp collisions. 

For energy loss we seek a form that is consistent, on one hand, with Ai? cx L as 
suggested in Ref. [97] for ID expansion, and on the other hand, with (dE/dL) cx E 
in Ref. [145 for 6 < i? < 12 GeV. A reasonable approximation of the differential 
energy loss is then 

AE 

— = PAL, (50) 

which translates to our variables as 

k-q = k(3t. (51) 

at small t with (3 being an adjustable parameter. For larger t we exponentiate the 
above to 

q = ke-^' , (52) 
and let G{q, k, t) take the simple form 

G{q,k,t) = q6{q~ke-^^^) . (53) 
The 5 function can be broadened to account for fluctuations, but we shall take 



( 53 1 as an adequate approximation of the complicated processes involved in the 
parton-medium interaction, the justification of which rests ultimately on how well 
the calculated result can agree with the data for 2 < < 12 GeV/c. 



Using ([531 in (|49| yields 

1 r'' 

which shows explicitly how fi{k) is transformed to Fi{q) by the nuclear effect 
parametrized by while fi{k) itself contains the hidden modifications due to 
such effects as iVcoii dependence and shadowing.!^ Since fi{k) is dNf^'^^ /kdkdy\y=Q, 



( 54 1 becomes in the limit of L ^ the invariant distribution for hard parton i pro- 



duction in pp collision. In heavy-ion collisions Eq. ( 11 ), written for centrality ~ 0%, 
can now be improved to the form 

S{qi) = j jF,{q)Si{q,/q) (55) 
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for any centrality with L being the geometrical path length. We can now write the 



TS recombination term in (48 1 as 



(56) 



where i?,^ is the RF given in ( 16 1 



Eq. (48) is now totally specified, relating the observable pion spectrum to the 



nuclear parameter /3L through Fi{q) given in (54 1. We emphasize that because all 



three components of TT+TS+SS recombination are included in ( 48 1 it can describe 



the pion distribution at all pt using only one parameter (3L for each centrality. Fig. 
11(a) shows the fits of the 7r° distributions for nine bins of centrality,'^^' with a 
different value of /3L used for each centrality.li^ Those values of f3L are shown by 
the nine points in Fig. 11(b) (ignoring the line for the moment). Evidently, the 
agreement with data is excellent over such a wide range of px and centrality c. 
Apart from concluding that the RM works well, the jet-medium interaction has 
been effectively summarized by one phenomenological function (3L(c), which could 
not have been extracted from the data without a reliable way to relate energy loss 
to the suppression of pion production at all px through appropriate hadronization. 
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Fig. 11. (a) Pion spectrEpH for various centralities lowered by a factor of 0.2 for each step of 
increase of 10%. The solid lines are fits in the RM^^'* with one parameter for each centrality c. 
Those parameters are shown by dots in (b), which are fitted by the curve f(c) that is the average 
of the dynamical path length ^ over P(^,c) defined in Eq. l|58|. (c) § distributions for various c. 



Having obtained PL{c) we now go a step further to inquire what kind of vari- 
ation of the dynamical path length can the nuclear overlap generate for any given 
centrality. That is, for a fixed c, (3L{c) is the average of a variable ^ over a proba- 
bility distribution P(^, c) that describes the likelihood that a particular trajectory 
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occurs at centrality c with an effective energy loss such that 

PL{c)^ I dCCP{C,c). (57) 

Thus ^ plays the role of /3L except that it can vary from to a maximum for every 
fixed initial elliptic geometry depending on the initial point and orientation of the 
hard parton. So ^ is the dynamical path length of a trajectory whose average ^(c) 
is f3L{c). In Ref. [144J P{C, c) is chosen to have the form 

^'(C,c) = ^f^(Co-0"^ (58) 

where N normalizes the total probability to 1. The two parameters Sixid a are 
adjusted to fit the nine points in Fig. 11(b), with the result shown by the curve 
that renders an excellent fit for 

^0 = 5.42, a = 15.2 . (59) 



Thus ( 58 1 is a very efficient way to describe the energy loss effect for all centrali- 
ties. The shapes of P{£,, c) are shown in Fig. 11(c) for 6 values of c, exhibiting the 
expected peak that decreases with increasing c or shrinking ellipse. In view of the 
difficulty of deriving f3L{c) from first principles, let alone P(^,c), it is very con- 
venient to have the path-dependent quenching effect be represented by the simple 



description in ( 58 1 and ( 59 1 



The single-pion distribution at midrapidity for centrality c can now be expressed 

as 

fd^PiLc)'-^, (60) 
pap J pap 



where (48l is to be used for dNT^{c,S,)/pdp, provided that Fi(q) in it is not as given in 



(54l, but with PL replaced by 5, and, of course, fi{q) scaled by A^cou(c). It should 



be clear that (48 1 has basically two parts: Fi{q) that describes the hard parton 
part and the rest that describes the hadronization part. We have in this subsection 
incorporated the centrality-dependent energy- loss effect on Fi(q) by the use of just 



two dimensionless parameters given in ( 59 ) . Having successfully formulated the 



treatment of the single-hadron distribution, we are now ready to proceed to the 
study of dihadron correlation in near- and away-side jets. 



6.2. Near-side and away-side yields per trigger 

For dihadron correlation many momentum vectors of partons and hadrons are in- 
volved. To depict clearly their relationships with one another, we show in Fig. 12 a 
pictorial representation of all of them. The near side is on the right and the away 
side on the left. The vectors k, q, and pt are, respectively, the momenta of the ini- 
tiating hard parton, of the same parton as it leaves the medium, and of the trigger 
hadron. The associated hadron on the same side is labeled by pa- On the away side 
the corresponding momenta are k' , q' , and pb, there being no trigger on that side. 
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If there are no transverse momenta in the beam partons, k and k' should be equal 
and opposite in every event; however, their averages over triggered events may not 
be the same, so it is better to label them distinctly from the start. The momenta of 
interest in the following are such that the hadronic px are in the range 2 < < 10 
GeV/c but with special emphasis on pT < 6 GeV/c. 




Fig. 12. A sketch of momentum vectors of partons (in red) and hadrons (in blue) with near side 
being on the right and away side on the left. 



For two pions in the same jet, neither of which are in the ridge, we can leave 
bution, and write as a g' 

PaY ^ J q I 



out the TT contribution, and write as a generalization of (48 1 

PtPadptdpa {PtPa 



TS(<z,p*) + ^AM- 



Pt 



TS(g -pt,pa) 



PtPa „^ 



ft Pa 

q ' q 



(61) 



where D2{zi, Z2) is the dihadron fragmentation function.^^ For notational brevity 
(61 1 is for a fixed ^, the averaging over P(^, c) being a process that can be applied 
when c is fixed. For two pions on opposite sides the recoil parton must be considered 
explicitly, so ( 49 1 should be generalized to 

r« (3dt 



Fliq^q'A)^ 



dkkf,{k)G{q,k,t)G{q',k',^-t) 
dkkf,{k)qq'5{qq' -kk'e-^). 



(62) 



with k' = k. Thus the dipion distribution for pf, on the side away from pt is 



dN^AO 
PtPbdptdpb 



—y 



dqdq' , , 
-yFM^q.i) 



TS{q,pt) + ^D: 



T5{q',pb) + ^D:, 



(63) 



The S function in ( 62 1 restricts the integration of q' in ( 63 ) to the range from qe 
to qe^ that correspond to the hard scattering point being on the near-side boundary 
to the far-side boundary. 
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The yield per trigger at fixed centrality can now be obtained for the near side 

as 

Y^^ (PuPa,c)- jd^p^^^c)dNAO/Ptdpt ' 

and similarly for Y^^'^^ {pt,pb, c) with dNT^jriO /PtPbdptdpb being used in the numer- 
ator. Since every factor is already specified, it remains only for the computation 
to be carried out. The results are shown in Fig. 13(a) for near side and (b) for 
away side. The three sheets for c — 0.05, 0.35, and 0.86 are separated by a factor of 
10 between sheets for clarity's sake. Superficially, the two figures may look similar 
in a vertical scale that spans over 4 orders of magnitude, but there are significant 
differences that only by closer examination can one learn from them the nature of 
the medium effects on the jets. 




Fig. 13. (a) Near-side jet yield per trigger vs trigger momentum pt and associated particle mo- 
mentum Pa- (b) Away-side jet yield per trigger vs pt and p;,- 



6.3. Medium effects on dijets 

Let us first cut the two figures in Fig. 13 by three fixed-pt planes at pt = 4, 6, and 8 
GeV/c, and show the results in Fig. 14(a) and (b) for c = 0.05 and 0.35. Note (i) the 
dependence of Y^^^'^{pa) on pa is more sensitive to pt than that of the dependence 
of Y^^^y{pb) on pb] (ii) the increase of the yield with pt is more pronounced for the 
near side than for the away side; and (iii) Y^^^" has negligible dependence on c, but 
y^way jjicj-eases by roughly a factor of 2 when c changes from 0.05 to 0.35. Let us 
discuss these three features separately. 

On item (i) one can determine the near-side average inverse slope Tq of the 
approximate exponential behavior in the range 2 < < 4 GeV/c for c — 0.05, 
for which there are data. The result is shown by the line in Fig. 15(a), and agrees 
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Fig. 14. (a) Left panel: near-side jet yield per trigger vs pa for fixed pt and c. Data points are 
from Ref. [6.11]. (b) Right panel: the same for away side vs p^,. 



well with the dataEin2lIIlZl Evidently, the spectrum of the associated particles in 
a jet becomes harder as the trigger momentum increases, as one would expect. It 
should be remarked that the line in Fig. 15(a) is for pions, while the data are for all 
charged particles. However, since pions dominate in jet peaks (unlike the ridges), 
the comparison is not unreasonable. The inverse slope T;, for the away-side jet at 
c = 0.05 for the same ranges of pt and pb is shown in Fig. 15(b), exhibiting a lower 
value compared to Ta and a mild decrease with pt. There exist no suitable data for 
comparison. The strong difference between Ta and has good reasons, as will be 
discussed below. Two data points from Ref. |148j are included in Fig. 14(a), the 
details of which are discussed in Ref. [144J. They lend support to the theoretical 
curve in both magnitude and pa dependences. 




0.4 



456789 456789 
(GeV/c) Pj (GeV/c) 

Fig. 15. Inverse slopes of associated particles on the (a) near side and (b) away side. Data in (a) 
are from Ref. [101] [147]. 



Note (ii) about the pt dependence of the yield is related to note (i) about the 
Pt dependencies of the shapes in pa and p^. First of all, Y^^^'^ is larger than Y^^^'^ 
in magnitude, meaning that there is more suppression on the away side than on 
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the near side. To quantify that interpretation let ^ be fixed at 2.9 corresponding to 
the maximum probability for c = 0.05 so that the suppression effect is not partially 
hidden by averaging over ^, which can vary from to ^q. Let the suppression factor 
for the near side be defined by 

-/3t\ 



r„car(pT) = {e-nPT = (9/fc> 



/Pt 7 



(65) 



where the average is performed over dN^^/pTdpT given in (48 1, and the last ex- 
pression follows from (53 1. r„car(PT) gives a measure of the fraction of momentum 



retained after energy loss reduces k to q. For the away side the suppression factor 
is defined by 



f{Pt,Pb) 



(66) 



away V/^t I ; — \i- IPt.Pb 

where — ^ + I3t is equivalent to —P{L — t) if ^ is denoted by f3L for a fixed medium 
length L, in which t is the portion from the hard-scattering point to the near-side 
surface, and L — t is the distance to the away-side surface. The average in (66) is 



done using the dihadron distribution given in ( 63 1 
1 

I- 

CL 

SO. 8 
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Fig. 16. Suppression factors for the (a) near side and (b) away side. 



The calculated results for rncar(PT) and raway(pt,P6) are shown in Fig. 16(a) 
and (b). They clearly indicate that rncar(PT) is much larger than raway(pt7Pti) with 
the implication that there is less suppression on the near side than on the away 
side. The physics of the phenomenon is clear: at large px the point of creation of 
hard parton with large k is predominantly close to the near-side surface in order 
to minimize energy loss with the consequence that the distance to the away-side 
surface is longer, thus more suppression for pb on that side. Since rncar(PT) saturates 
at around 0.85, only 15% of the parton energy is lost to the medium on the near 
side. The corresponding {f3t) is less than 0.2, so {(it)/S^ — {t) / L « 0.065, meaning 
that the hard partons are created within a layer of thickness ^ 13% of L from the 
surface. That is a quantitative description of trigger bias. On the other hand, the 
behavior of raway(PtiPfc) reveals the opposite: at fixed pf, it decreases with increasing 
Pt, implying more suppression as the hard-scattering point is pulled closer to the 
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near-side surface. That is antitrigger bias. At fixed pt, raway(pt,Pb) increases with 
Pb, since higher pb demands higher q' , which can be satisfied only if {L — t) is 
reduced or (fc') increased, actually both. Lowering (L — t) is, of course, a way to 
reduce energy loss by pulling the scattering point closer to the away-side surface. 
But the increase of {k') with pb is another aspect of antitrigger bias, whose details 
are described in Ref. |144j . Event-by-event momentum conservation requires k' = k; 
however, when averaged over all events, (fc') depends on pb, while (fc) does not. In 
general, (fc') is far greater than (k) because any trigger favors shorter (t), and any 
finite pb pushes up {k')/{k) though not {k' /k). Because of the difficulty of producing 
an associated particle on the away side relative to one on the near side, Tb is lower 
and decreases with increasing pt in Fig. 15(b), while Ta is higher and increases with 
Pt in Fig. 15(a). 

Now on note (iii) about centrality dependence Fig. 14(b) shows that the depen- 
dencies of Y^^'^y{pt,pb,c) on Pt and pb are essentially the same whether c = 0.05 
or 0.35, except that the magnitude of the yield is increased due to the reduction 
of path length at higher c. For the near side there is essentially no dependence of 
Y^^'^'^{pt,Pa, c) on c. That is shown more explicitly in Fig. 17, where the lower three 
lines are for pt — A GeV/c and the upper line is for pt — Q GeV/c. The solid lines 
are for integrated yields with 2 < < 4 GeV/c. The near independence on c is 
a manifestation of the trigger bias, since the hard-parton production point, being 
restricted to a layer roughly 13% of L just inside the near-side surface, is insensi- 
tive to how large the main body of the medium is. Actually, the decrease of the 
TS component with c balances the increase of the SS component with c so that the 
net yield being their sum is approximately constant in c. The data points in Fig. 
17 support the calculated result in both the magnitude and the c dependence of 
l^^^''(pt, c) when pa is integrated from 2 to 4 GeV/c^^ 
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Fig. 17. Near-side yield per trigger vs centrality c. The solid lines are for integrated yields with 
2 < po < 4 GeV/c. The data points are from Ref. |101| . 
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In summary the discussion above gives quantitative demonstration that the 
trigger bias is the preference for the hard-scattering point to be close to the near- 
side surface and that the antitrigger bias is the consequence: (fc') is much larger 
than (fc), and even larger than pt or pi,. Those are the properties of the events 
selected by a trigger at pt with an associated particle on either the near or away 
side. 



6.4. Symmetric dijets and tangential jets 

We now follow the summary comment at the end of Sec. 6.3 with the question on 
what if we select events with symmetric dijets where pt = Pb- Instead of studying 
the properties of a third particle in association with the two trigger particles (for 
which we need trihadron correlation function), we can nevertheless learn a great 
deal from examining closely various calculable quantities in the dihadron correlation 
problem. Let p be the momentum of the symmetric dijets, p = Pt = Pb- One can 
calculate {k'){p,c) and {q'){p,c) for various centralities, even though they are not 
directly measurable. The result is that they both increase almost linearly with p 
and that there is essentially no dependence on c. Their ratio {q') / {k') is therefore 
approximately constant in p with a value of about 0.8 as shown in Fig. 18(a). 
That behavior is similar to the property of rnoar(p), shown in Fig. 16(a), which, 
according to (65 1, is also (q/k). It is important to bear in mind that {q') and (fc') are 
averages over the two-particle distribution dNT^Tr{c) / PtPbdptdpb with pt — Pb, given 
in (63) for fixed ^, followed by averaging over P(^,c), whereas (q), (fc) and {q/k) 
are averages over the single-particle distribution dN^^/pdp, as noted after (65 1. The 
near-side averages know nothing about the away-side analysis, so r,ioar(p) describes 
only the suppression associated with trigger bias. The fact that the suppression on 
the away side is about the same as on the near side, when the average (q'/fc') is 
over symmetric dijets and (q/k) is averaged over near-side jet, and that it is true 
for any centrality, has only one important implication: the dijets are created very 
near the surface on both sides. It means that the symmetric dijets are dominated 
by tangential jets, which behave similarly at all centralities. That is a striking 
conclusion. Hard scattering can, of course, occur anywhere in the overlap region. 
But those partons from the interior lose most of their momenta on the way out. 
Those that are created on the near side would not give rise to symmetric dijets. 
Thus only those created along the rim and directed tangentially to the surface can 
lead to events with pt = Pb- 

There is some experimental evidence to support this finding. In 2jet + 1 correla- 
tion studiecPal it is found that the third particle does not show any ridge structure 
and that the centrality dependence goes as iVpart- The latter means that the dijets 
are created near the surface, while the former means that the dijets are tangential 
because we have already seen in Sec. 4.4 that ridge production depends on the 
local flow direction to match the trigger direction due to the correlation function in 
(27 1. The flow direction near the surface is normal to the surface and is therefore 
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P (GeV/c) (GeV/c) 

Fig. 18. (a) The ratio {q')/{k') at the symmetry point p = pt = Pb for four values of centrahty. (b) 
Distribution of associated pion (p;,) in the away-side jet for six values of pion trigger momentum 
(pt) in GeV/c for ^ fixed at 2.9. 

normal to tangential jets. In the experimental analysi^^the two jets have pt cuts 
at > 5 GeV/c and px^ > 4 GeV/c, so they are not exactly symmetrical. If ridge 
formation from those dijets is to be discovered in the future, it would pose a serious 
challenge to the treatment of ridgeology in Sec. 4. But so far the data are in accord 
with the results both here and on ridges. 

6.5. Unsymmetric dijets and tomography 

The conclusion in the last section on the dominance in symmetric dijets by tangen- 
tial jets leads one to consider the only option left for probing the dense medium 
by parton jets apart from using direct 7, and that is the study of away-side jets in 
unsymmetric dijets. If hard or semihard partons created near the surface are most 
responsible for the near-side jet, then the away-side jet should experience fully the 
effect of the medium. That is indeed the case when one calculates {q') and (fc') 
for various values of pt and pb at fixed ^. It is found in Ref. |144| that for ^ = 2.9, 
{q') /{k') is in the range of 0.15 to 0.3 for 4 < pt < 10 GeV/c and 2 < pfc < 6 GeV/c, 
and gives a quantitative measure of the degree of energy loss. Unfortunately, that 
cannot be checked experimentally, since only centrality can be selected in realistic 
collisions, not the dynamical path length ^. 

At fixed centrality the value of ^ can fluctuate over a wide range, as can be 
seen in Fig. 11(c). For c = 0.05 the average ^(c) is 2.9, but the characteristic of the 
away-side jet is dominated by the lower ^ portion of the range. The consequence 
is that Y™^y{pt,ph,c) does not depend sensitively on pt, as Fig. 14(b) has already 
indicated for any fixed c. That is in sharp contrast from the case where ^ is fixed 
at 2.9, and Y™'^y{pt,piy,^) decreases by an order of magnitude as pt is decreased 
from 9 to 4 GeV/c at any fixed pt, as shown in Fig. 18(b). That diminishing yield is 
because (fc) is lower at lower pt and the energy loss by the recoil parton traversing 
the thick medium results in reduced probability of producing a pion at fixed pi,. 
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That is not the case in Fig. 14(b). At fixed c, the decrease of pt does not lead to 
significantly lower Y^^'^^{pt,Pb,c) because the hard-scattering point is already in a 
region of lower ^ to minimize energy loss. Lowering pt increases the yield at fixed 
Pb due to softer parton, but the number of triggers is also higher, so the yield per 
trigger remain nearly unchanged. In other words, allowing ^ to be small even at 
small c removes energy loss as a decisive factor in the problem. That is what makes 
dijet tomography ineffective as a probe to learn about a medium that has no fixed 
thickness. 

Further insight can be gained by studying {(3t') at fixed c, where t' denotes 
the distance from the hard-scattering point to the away-side surface. Fig. 19(a) 
shows {f3t') vs Pb for various pt and c^^^ A general impression from that figure is 
that iPf) is low, less than approximately 0.4. That is much lower than PL = 2.9, 
which is the average dynamical path length determined from fitting the single- 
pion distribution from c = 0.05 [see Fig. 11(b)]. The height of that distribution 
at large pt [see Fig. 11(a)] is what renders Raa ~ 0.2, a number that pQCD 
calculations aim to obtain. The fit in Ref. |144j is achieved by setting /3L = 2.9 
in order to obtain the correct normalization for (INt^ /ptdpt at large p^, for which 
the contributions from all partons, near or far, hard or semihard, are counted. The 
result that {(3t') <C (iL shown in Fig. 19(a) is therefore indicative of the fact that 
conditional probability with pt and pb fixed is highly restricted compared to the 
inclusive probability. With {t') being much less than L at fixed c, one is led to 
conclude that the unsymmetric dijets are also produced near the surface and are 
essentially tangential, as with symmetric dijets. Thus the medium interior is not 
probed. This conclusion is distinctively different from the case where ^ is fixed. Fig. 
16(b) has shown that the suppression on the away side can be large (hence Faway 
small) for unsymmetric dijets at ^ = 2.9. In that case {(it') would not be small, as 
we now have for any fixed c. 
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Fig. 19. (a) Left panel: average dynamical path length of recoil parton directed toward the away 
side, (b) Right panel: Ratio of the average parton momenta of recoil parton at the surface to that 
at the creation point. 
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Experimental data-'^'*^ on slightly asymmetric dijets at high pT show no essential 
difference in the structures of the two jets and arc in support of the findings in Rcf. 

Another way to come to the same conclusion is to study the energy loss of 
the away-side jet. In Fig. 19(b) is shown {q')/{k') for c = 0.05 and 0.35.IS1 It is 
the fraction of average momentum retained by the hard or semihard parton after 
traversing the medium on the away side. Being around 0.8 implies that only about 
20% of the parton momentum is lost to the medium, not much more than the 
fractional energy loss on the near side. At fixed pt Fig. 19(b) shows an increase 
of {q')/{k') with increasing pti because the hard scattering point is pulled more to 
the away side, but it shows a decrease of the ratio at increasing pt since the point 
is then pulled to the near side. This push-and-puU effect of pt and is clearly 
what one expects in the oppositely-directed jets when the path lengths on the 
two sides are comparable. At the symmetry point pt = Pb — ^ GeV/c, {q')/{k') 
for c = 0.05 is only slightly lower than its value for c = 0.35, implying strongly 
that the fractional energy loss on both sides remains about the same regardless of 
centrality. That can only mean that the hard partons are created near the surface 
and directed tangentially. Making pt ^ Pb in unsymmetric dijets does not change 
{q')/{k') drastically. Thus so long as the medium thickness cannot be controlled, 
there seems to be no useful tomography that can be done with parton-initiated 
dijets. If the away-side jets are dominated by those created near the away-side 
surface (such as the tangential jets) , then the events triggered by direct 7 on the near 
side are also likely to be dominated by those where the hard scattering takes place 
near the away-side surface (not just tangential) so long as an associated particle 
with significant pb is required on the away side. The 7r-triggered and 7-triggered 
distributions, Iaa, of the associated particles should be roughly the same. There is 
some experimental evidence for that similarity.'^^ 

It is important to note that the above comment is for tomography only, i.e., 
medium effect on jets. The jet effect on medium, such as Mach cone, is a difference 
matter that depends on different physics and may well reveal properties of the 
medium interior. 

7. Conclusion 

In this review many problems in heavy-ion collisions have been examined in many 
parts of the phase space. The good agreement between theoretical calculation and 
experimental data in almost all cases cannot but affirm that the successes of the 
theoretical approach adopted cannot be all fortuitous and that the interpretations 
given to the physical origins of the measured phenomena are not without some de- 
gree of realism. In some cases there are no apparent alternative schemes to explain 
the data. Of course, a phenomenological model is not a theory based on first prin- 
ciples, but its usefulness should not be overlooked when its scope is out of reach by 
any theory commonly accepted. 
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In the first part of this review, mainly in Sees. 2 and 3, we have summarized 
the essential basis of the recombination model, answering some questions raised 

by critics, and pinning down some of the parameters, while pushing the frontier 
to kinematical regions no other models have attempted, and establishing parton 
recombination as a universal hadronization mechanism. Then in the second and 
main part recombination is used in an essential way to relate what are observed 
about jets and ridges to two complementary aspects about the dense medium: the 
effects of hard or semihard partons on the medium and the converse. Those aspects 
of the medium effects cannot be made empirically relevant unless there is a reliable 
description of hadronization for all pT where correlation data exist. 

Some of the new insights gained that are of particular current interest are: 

(a) correlation between trigger and ridge (Sec. 4.4), 

(b) independence of 4> anisotropy on fast thermalization (Sec. 5.1), 

(c) recombination does not guarantee quark number scaling (Sec. 5.3), 

(d) different properties of the near- and away-side jets (Sec. 6.3), 

(e) dominance of tangential jets in symmetric and unsymmetric dijets (Sees. 6.4 

and 6.5). 

If there is to be one unifying conclusion to be made as a result of these findings, 
it is that most observables on jets and ridges are due to hard or semihard partons 

created near the surface. Partons created in the interior of the medium that lead 
to dijets are not able to compete effectively with those that have shorter distance 
to traverse. So long as the observables allow the parton creation points to include 
regions that offer the partons a choice of paths of least resistance, they will take it 
and dominate. That is the reality faced by experiments that can only fix centrality, 
not medium thickness. With that recognition the efficacy of jet tomography is called 
into question. Unlike X-ray scanning of organic or inorganic substances, there is no 
control of the sources of the hard partons, so the necessary averaging process under- 
weighs the contribution from the region of the medium that one wants to learn most 
about. 

The above comment refers to parton-initiated dijets. Of course, for single hadron 
at large px the nuclear modification factor Raa has long been used as a measure of 
energy loss in dense medium in experiments and in theory. As soon as a condition 
is imposed on the detector of another particle on the away side, the region of the 
system probed in changed. If that hadron's pT is low, it may be in the double- 
humped shoulder region, which does provide some information on what the effect 
of the away-side jet is on the medium. However, to learn about the effect of the 
main body of the medium on the jet from dihadron correlation is more difficult. 

Ridgeology addresses a different set of problems, quite distinct from dijets. 
Ridges are stimulated by jets, but are not a part of the jets that are character- 
ized by the participation of shower partons. Considerable attention has recently 
been drawn to the study of ridges. At this point there is no consensus in their 
theoretical interpretation. The connection between the ridges found in triggered 
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event£^ili^ and those found in autocorrelation without trigger^MHU is in our 
opinion tenuous, until the pr dependence of the ridge is clarified. Minijet^^^ peak- 
ing at Ut ^ 2.8 correspond to (px) ~ 1.2 GeV/c, which is much lower than the 
Pt range of the triggers used in ridge analysis. The ridge found in autocorrelation 
has no pT cut, while those found in triggered events have p^^°'^ > 2 GeV/c.'^^ 
Nevertheless, there may exist a connection between autocorrelation and the ridges 
without trigger discussed in Sec. 4.2. 

Another phenomenon of some current interest is the observation of extended 
ridge at large Arj)^^ Before one concludes that such a long-range correlation can 
only arise from the mechanism of strings (or color flux tubes) being stretched be- 
tween forward- and backward-going quarks, it seems prudent to allow firstly the 
possibility of other types of early-time dynamics, and secondly a broader view of 
the problem of particle production at large \r]\. In particular, one should consider 
the role played by the hard parton that leads to both the trigger and the ridge. 
Moreover, one should consider the issue of large p/ir ratio at t] = 3.2 (see Sec. 3.4) 
as a part of the solution of the bigger problem. If the observed large p/w ratio is 
a feature of the final-state interaction (FSI), one should not regard the ridge phe- 
nomenon as a manifestation of the initial-state interaction (ISI) only without taking 
into consideration also the effects of FSI. As we have noted at a number of places 
throughout this review, hadronization is an important link between the observables 
and partonic dynamics. The characteristic of ridge formation in azimuthal correla- 
tion discussed in Sec. 4.4 is a good example of the interplay between ISI and FSI. 
It would be surprising if the same does not hold true for the correlation in rapidity. 

The considerations given in this review to jets and ridges may not be directly 
relevant at LHC where jets are copiously produced — unless px is extremely high. 
At RHIC the background to a rare jet is thermal, but at LHC the background to 
a high-pT jet, say at px ^ 100 GeV/c, includes many other lower-p^ jets. The 
admixture of thermal and shower partons in the background introduces new com- 
plications to the notion of enhanced thermal partons, and renders what is simply 
conceived for RHIC inadequate for LHC. That new energy frontier will indeed open 
up a wide new horizon. 
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